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We constrain cosmological models where the primordial perturbations have both an adiabatic and 
a (possibly correlated) cold dark matter (CDM) or baryon isocurvature component. We use both a 
phenomenological approach, where the power spectra of primordial perturbations are parametrized 
with amplitudes and spectral indices, and a slow-roll two-field inflation approach where slow-roll 
parameters are used as primary parameters determining the spectral indices and the tensor-to-scalar 
ratio. In the phenomenological case, with cosmic microwave background (CMB) data, the upper 
hmit to the CDM isocurvature fraction a is 6.4% at fc = 0.002 Mpc~^ and 15.4% at fc = 0.01 Mpc"^ 
At smaller scales (larger k) larger isocurvature fractions are allowed, and therefore large values of 
the isocurvature spectral index, riiao ~ 2, are formally favored. The median 95% range for the 
non-adiabatic contribution to the CMB temperature variance is —0.030 < ar < 0.049. Including 
the supernova (or large-scale structure, LSS) data, these limits become: a < 7.0%, 13.7%, and 
-0.048 < QfT < 0.042 (or a < 10.2%, 16.0%, and -0.071 < ar < 0.024). The CMB constraint on 
the tensor-to-scalar ratio, r < 0.26 at fc = 0.01 Mpc"'^, is not affected by the nonadiabatic modes. In 
the slow-roll two-field infiation approach, the spectral indices are constrained close to 1. This leads 
to tighter limits on the isocurvature fraction, with the CMB data a < 2.6% at fc = 0.01 Mpc~^, but 
since the non-adiabatic contribution to the CMB temperature variance comes mostly from larger 
scales its median 95% range is not much affected, —0.058 < qt < 0.045. Including supernova 
(or LSS) data, these limits become: a < 3.2% and -0.056 < ar < 0.030 (or a < 3.4% and 
—0.063 < ar < —0.008). When all spectral indices are close to each other the isocurvature fraction 
is somewhat degenerate with the tensor-to-scalar ratio. In addition to the generally correlated 
models, we study also special cases where the adiabatic and isocurvature modes are uncorrelated 
or fully (anti)correlated. We calculate Bayesian evidences (model probabilities) in 21 different cases 
for our nonadiabatic models and for the corresponding adiabatic models, and find that in all cases 
the current data support the pure adiabatic model. 

PACS numbers: 98.80.-k, 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 



The nature of primordial perturbations is the key ob- 
servation window to the physics of the very early Uni- 
verse and the very-high-energy physics that is beyond the 
reach of particle accelerators. The simplest model for the 
generation of primordial perturbations, single-field slow- 
roll inflation, predicts adiabatic Gaussian nearly scale- 
invariant primordial scalar perturbations, and tensor per- 
turbations at a level that may or may not be observable, 
depending on the detailed model. Current cosmological 
data are consistent with this kind of perturbations, with 
so far no evidence for primordial tensor perturbations. 
There is some evidence of a small deviation from scale 
invariance. More complicated models, like multi-field in- 
flation, may give rise to deviations from adiabaticity or 
Gaussianity that improved observations may be able to 
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detect. 

In this paper we study what information currently 
available data on the cosmic microwave background [IH3] 
give on possible deviations from adiabaticity, and check 
how the results change if we add supernova [4] or large- 
scale structure [5] (that includes baryon acoustic oscil- 
lations, BAG [S]) data into the analysis. While adia- 
batic primordial perturbations are fully determined by 
one physical perturbation quantity, which we may take to 
be the curvature perturbation in the comoving gauge, TZ, 
deviations from adiabaticity represent additional, isocur- 
vature, degrees of freedom that can be given in terms 
of entropy perturbations, which represent how perturba- 
tions of different particle species differ from each other. 
For example, the primordial cold dark matter (GDM) 
entropy perturbation can be written as 



Sc-f — 



(1) 



where ric and are the number densities of CDM parti- 
cles and photons at early time, deep in the radiation dom- 
inated era. This kind of perturbations do not (initially) 
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perturb the spatial curvature of comoving slices {TZ = 0), 
hence the name "isocurvature perturbations". There is 
so far no evidence for an isocurvature perturbation com- 
ponent, and observations require it to be subdominant 
compared to the adiabatic component [THU] • Even if no 
isocurvature component is found, improving upper limits 
to it will help to constrain the particle theory related to 
the origin of perturbations [20^. 

Out of four possible isocurvature modes [HI [55] , we fo- 
cus here on cold dark matter (CDM) and baryon isocur- 
vature perturbations, which are observationally indis- 
tinguishable from each other (see however [23] for the 
prospects of using the Hydrogen 21 cm Hues), and do 
not consider the less well motivated neutrino isocurva- 
ture perturbations [21] • In our earlier work [5J (TT] |5S] we 
assumed that tensor perturbations were negligible. Now 
we allow for the presence of tensor perturbations. We 
assume a power-law power spectrum for the primordial 
curvature and entropy perturbations and allow for corre- 
lations between them. We also study special cases with 
no correlation or with ±100% correlation. 

Assuming spatially flat geometry of the Universe, we 
perform full parameter scans of this mixed (adiabatic and 
isocurvature) model, as well as the standard pure adia- 
batic ACDM model for comparison. We present posterior 
probability densities of the standard cosmological param- 
eters and the extra isocurvature parameters, and report 
the Bayesian evidences for these models calculated with 
the MultiNest nested sampling package p6], using the 
CMB data alone or CMB and supernova data or CMB 
and large-scale structure data. 

We use two different approaches: 1) A phenomenolog- 
ical approach, where we make no reference to the origin 
of the primordial perturbations, and just try to deter- 
mine or constrain their amplitudes and spectral indices 
from the data. 2) A slow-roll two-field inflation approach, 
where we assume the perturbations were generated by 
quantum fluctuations during two-field inflation, and the 
spectral indices are determined by the slow-roll param- 
eters at the time the cosmological scales exited horizon 
during inflation. 



II. CURVATURE AND ENTROPY 
(ISOCURVATURE) PERTURBATIONS 

In this section we deflne the curvature and isocurvature 
perturbations, as well as our notation and sign conven- 
tion for them. 

In general, the metric perturbed by scalar perturba- 
tions is given by the line element 

ds'^ = - (1 + 2(/))dr2 + 2d^B drdx' 



(1 - 2ijj)S,j + 2didjE 



dx^dx^ 



} 



(2) 



gauge- invariant Bardeen potentials [57] : 

$ = 4> + n{B-E') + {B- E')' , (3) 
^ = ij~HiB~E'), (4) 

where H = a' /a is the conformal Hubble parameter, a 
is the scale factor of the Universe, and prime denotes 
a derivative with respect to conformal time r. In the 
longitudinal (Newtonian) gauge, where B = E = 0, the 
Bardeen potentials match the Newtonian metric pertur- 
bations. 

The comoving gauge is deflned by requiring both co- 
moving slicing, B = 0, and comoving threading, u = 0, 
where v is the fluid velocity potential. We denote the co- 
moving gauge by superscript C, and define the comoving 
curvature perturbation TZ by 



n 



— j.c 



(5) 



Here we used the gauge transformation properties of ^ in 
order to relate ^/j*^ to the perturbation quantities in any 
other gauge G. 

Our sign convention for TZ is such that a positive TZ 
corresponds to positive curvature of the 3-dimensional 
T = constant slice. If we go to a spatially flat gauge 
(superscript F) where "ip^ — E^ — 0, we flnd 



TZ = 



3(1 + w) 



5p^ 
P 



2 / fc 

3 VH 



(6) 



where p is the background (average, unperturbed) energy 
density, and w — p/ p with p the background pressure. A 
positive TZ corresponds to an overdensity, 5p^ > 0, on 
large scales k/Ti <C 1. 

Finally, the definition ([5| can be expressed in various 
forms in terms of the gauge invariant Bardeen potentials: 



TZ 



3(1 + w) 



(■H-^*' + $) 



(7) 
(8) 



It is important to note that our sign convention is op- 
posite to, e.g., [55H5n] . and whether appear in the 
time or space part of the metric and what sign they 
have in Eq. ^ varies a lot in the literature. Further- 
more, often the definitions of ^ and $ are swapped and 
also their sign convention varies. Our sign convention 
matches with our previous works [Q] HT] [25] and, e.g., 
with 19J. (Ref. [2D] uses our sign convention for TZ, but 
has ^ and $ swapped.) 

We define the total entropy perturbation by 



S = V. 



5p Sp 
p' p' 



(9) 



From the above metric perturbations one can construct 



where p' and p' are the conformal time derivatives of 
the background pressure and energy density, respectively. 
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Further, we define an entropy perturbation between par- 
ticle species x and y by 

' 6px _ 5py^ 

P'x P'y 



(10) 



Using the continuity equation, this can be written as 
Sxy = + Wx) - 5yl{l + Wy), where 5 = 5pl p. This 

leads easily to Eq. ([T]) for CDM and photons. Both S 
and Sxy are gauge invariants, and hence we did not need 
to specify the gauge above. 

It can be shown that on super-Hubble scales (fc/'H <C 
1) U-^n' oc S [lll2il3IH33]- Therefore, in the absence 
of entropy perturbation, the comoving curvature pertur- 
bation remains constant on super-Hubble scales. In the 
contrary, 7i~^S' cx S, and hence the total entropy per- 
turbation evolves with time even on super-Hubble scales. 
Moreover, these differential equations imply that (part 
of) S can be "converted" into TZ, but TZ cannot be con- 
verted into S, and even if TZ and S were initially (e.g. at 
horizon exit during inflation) uncorrelated, the later evo- 
lution can lead to a correlation between them: there will 
be a part of TZ that has been "created" from the initial 
iS, and hence these are fully correlated or anti-correlated 

If we assume that the Universe consists of photons 
(7), neutrinos (i/), baryons (6), and cold dark matter 
(c) , and that the entropy perturbation between neutrinos 
and photons vanishes, then the total entropy perturba- 
tion during radiation domination can be written as 



S = 



(11) 



where Ueq is the scale factor at radiation-matter equality, 
fc Pc/iPc + Pb), and fb = Pb/{Pc + Pb)- 

In this paper, we perform all the isocurvature analysis 
by specifying a non-zero primordial Sc-f, i.e., a primordial 
cold dark matter isocurvature mode, while keeping all the 
other relative entropy perturbations S^y zero. From now 
on we shorten our notation, and denote 



S=S, 



C7 



(12) 

at primordial time trad, deep in the radiation dominated 
era, on super-Hubble scales. 

Finally, we comment once more on our sign conven- 
tion. The ordinary Sachs- Wolfe (SW) effect in the CMB 
temperature anisotropy, in the direction n at the sky, is 

ST 1 

— (n) = - - [7^(^rad,n) + 2/,5(t,ad,n)] . (13) 
1 5 

(Note that we have —1/57?. for the adiabatic part whereas 
equation (3) in [33 shows -fl/57?. due to their different 
sign convention for TZ.) Note that TZ and S appear with 
the same sign in this expression. Hence, a positive pri- 
mordial correlation, C^s > 0, leads to extra power at 
low multipoles in the CMB temperature angular power 
spectrum, Cf cx {\5T/T\^), compared to the uncorrelated 
case. Vice versa, a negative primordial correlation leads 
to a negative contribution to the SW effect, and hence 
reduces the angular power at low multipoles. 



III. MECHANISMS THAT PRODUCE 
ISOCURVATURE 

There are various mechanisms that can produce (cor- 
related) isocurvature and adiabatic perturbations. Com- 
pared to the single-field slow-roll inflation which can 
stimulate only the adiabatic perturbation mode, at least 
one extra degree of freedom is needed. Inflationary mod- 
els with multiple (scalar) fields are a natural extension. 
Production of isocurvature perturbations in this type 
of models has been studied, e.g., in Refs. [251 EOl ISH" 
[7T] . some references going to a non-linear level [72H74j . 
Many of these references derive so called consistency re- 
lations (there exists actually an infinite hierarchy of them 
[751 [75]), or analytical expressions for the spectral tilts 
and amplitudes of adiabatic and isocurvature compo- 
nents in specific inflationary set-ups. We will use the 
notation and the leading order results of [3D] for two-field 
slow-roll inflation. 

Various types of axions carry isocurvature perturba- 
tions [771190) ■ and it is possible to construct axion mod- 
els which lead to a large isocurvature spectral index, 
TT-iso ^ 2 — 4, [21]. (These values are formally preferred 
in our phenomenological approach.) Most axion-dilaton- 
moduli models give rise to isocurvature perturbations 
|92) . Moduli fields and isocurvature have been studied 
]. One solid candidate for producing correlated 



mixtures of adiabatic and isocurvature perturbations is 
the curvaton model [371 [25], or curvaton-type models, 
such as late-decaying scalar condensations [M]- Isocur- 
vature perturbations in these models have been discussed 
in Refs. jlOOHllS] . A generic feature of (minimally) su- 
persymmetric (standard) model is the existence of ex- 
tra fields that can carry isocurvature perturbations [TT^J- 
1122] . Refs. [I23lfl29] focus on the AfHeck-Dine mecha- 
nism which typically leads to isocurvature perturbations. 
In brane models isocurvature also arises naturally [1301 - 
I134j as well as in the ekpyrotic scenario [1351 1136] . 

Dark energy isocurvature perturbations have been 
studied in [137, . In some models with interaction between 
CDM and dark energy, an isocurvature mode is rapidly 
growing in some regions of the parameter space overtak- 
ing the usual adiabatic mode, and hence these regions 
can be ruled out |138H141] . Isocurvature in quintessence 
cosmologies has been studied in |142H147| . Large-scale 
magnetic fields can also induce a magnetic isocurvature 
component |148H150] . Finally, cosmic strings and other 
topological defects generically create isocurvature pertur- 
bations, but in this case they may not be "primordial" 
but instead continuously created during the evolution of 
Universe |151H155] . In order to match cosmic string mod- 
els with observations one needs another mechanism to 



produce the predominant adiabatic component: in |156j 
a model with a mixture of adiabatic perturbations from 
inflaton, and isocurvature perturbations from strings was 
observationally constrained. 

After introducing our model and its parametrization 
in the next section, we give a few specific examples in 
Sec. UVCl 
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IV. THE MODEL 

The general perturbation, here presented in Fourier 
space, is a superposition of adiabatic and isocurvature 
perturbations, 7?.(k) and S'(k). Its power spectrum V = 
Vn+Cizs +0371+^^3, where C represents the correlation, 
is defined by the expectation value 



[7^(k) + 5(k)]*[7^(k) + 5(k)]^ = (2^)3,5(3) (k-k)x 



fc3 



[Vn{k)+Cns{k)+Csnik)+rsik)] . (14) 



Note that we used a calligraphic letter V, while P is 
reserved for P{k) = ^V{k). 



Following m nil 1251 [TFT] we divide the comoving 
curvature perturbations into an uncorrelated part ( "ar" ) , 
and a part fully correlated with the entropy perturbation 
("as"), and assume power-law forms for the spectra: 



7'K(fc) -Par(fc)+7'as(/^) 



where 



and 



Par(fc) - v4'o(fc/fco)' 



(15) 

(16) 
(17) 
(18) 



Vs{k)^Bl{k/k^r--' . 
The covariance is given by 

Cnsik) = Csn{k) = ^oi?o (fc/fco)"™'"' , 
where 

ncor = {rias + n-iso)/2 . (19) 
In general, we denote ^-p^^^ki), A^- = 'Pa,s{ki), and 

For fitting the different models to the data, we choose 
our reference scales as 



fci = 0.002 Mpc" 

fco = 0.010 Mpc" 

fca = 0.050 Mpc" 

We further define 

l2 _ ^2 1 a2 , d2 „, _ d2 / 42 



(20) 



4^ = AU + Af, + Bt , a, = B/M^ (21) 



and 



7i = sign(A5jBi) 



(22) 



so that 

i2 _ n I , \\r a2 , /i2 



^ il^\j,\)iAt.,+Ai,) = (l-|7.|)(l-a.)A^ 



B? 



\l^\iA'„+Al 



H\{\-a,)A\ 



AsiBi = aco^Ai =Cns{ki)=Csn{ki), 



(23) 



where we defined the relative amplitude of the primordial 
correlation between the adiabatic and CDM isocurvature 
perturbations by 



sign(7i)^Q;i(l - Q;i)|7i 



(24) 



The total CMB temperature angular power spectrum can 
be written as 

Q - Al[{l - ao)(l - boDCT + (1 - ao)|7o|CT 
+ aoCf ° + acoroC'™'' + (1 - Q;o)roCj] 



Car I y^as < /^iso i i^cor > /^T" 
^ -t- -h -r 0£ -h , 



(25) 



where the Cg represent the different contributions to the 
angular power spectrum that would result from a corre- 
sponding primordial spectrum with a unit amplitude at 



the pivot scale k = ko (see 



The last term, Cj , is 



the possible contribution from primordial tensor pertur- 
bations, see the next subsection. 

From Eq. (11) we notice that in order to lead to 
the same observational effects we would need a pri- 
mordial baryon isocurvature perturbation that is Sty — 
{ fclfb)Sc-i = (j-^c/^b)S. So, if we find a constraint Be for 
Bi (in case of cold dark matter) , this can be converted to 
the constraint Bb — {ujc/^b)Bc for the baryon isocurva- 
ture. Therefore a constraint ac for ai can be converted 
to a constraint 



ab = 



1 



(26) 



For example, if ujb = 0.02 and oJc — 0.10, then ac — 0.10 
corresponds to at ~ 0.74. Unless otherwise stated, as 
in the following refer to the CDM isocurvature fraction. 
The constraints on the primordial baryon isocurvature 
fraction reported in this paper have been obtained from 
our runs with CDM isocurvature by the above mapping 
and by weighting the posterior by a Jacobian determinant 
of the mapping (in this case simply \dab/ dad) , i'^ order 
to simulate a situation where ab would have had a fiat 
prior 01171153]. 

Additional derived parameters can be defined. For ex- 
ample, a parameter n'^ represents the spectral index for 
adiabatic modes obtained expressing the adiabatic con- 
tribution as a single power law: 



d\nk 

(nar-l)(l-|70|)^""-' + (nas-l)|7|fc"" 



(l-|7o|)fc" 



|7o|fc'^ 



(27) 



where k = k/kg. 

Our pivot-scale free measure of the non-adiabaticity 
will be the total non-adiabatic contribution to the CMB 
temperature variance 

(((5r"°"-^'^)2) 



aT 



^^^^total from scalar pcrturbatioiis^2^ 

Efc2"(2^ + i)(cr + CT + cf° + Q") 



.(28) 
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Parameter 



Explanation 



range (min, max) 



Primary background parameters (common for both parametrizations) 



1000 

T 



physical baryon density; oj^ = /i Jlj, 
physical cold dark matter density; ujc = h^flc 
6 is the sound horizon angle; 6 — rs{zf) / Da{z,) 
optical depth to reionization 



(0.010, 0.050 
(0.02, 0.30 
(0.5, 2.2 
(0.02, 0.30 



Primary perturbation parameters in amplitude parametrization 



ln(10iM?] 

71 

I72 

ai 

Q-2 

( '■0 ) 
Asz 



Ai is the overall primordial perturbation power at = fci = 0.002 Mpc~^ 
A2 is the overall primordial perturbation power at k — k2 = 0.05 Mpc~^ 
primordial ratio of correlated adiabatic component to total adiabatic power at = fci 
primordial ratio of correlated adiabatic component to total adiabatic power a,t k — k2 
primordial isocurvature fraction at fc = fci; Alai = primord. isocurvature power at ki 
primordial isocurvature fraction at fc = fc2; ^202 = primord. isocurvature power at fc2 
primordial ratio of tensor and scalar power spectra at fc = fco = 0.01 Mpc~^ 
amplitude of the Sunyaev-Zel'dovich (SZ) template 



(1.0, 7.0 
(1.0, 7.0 
(-1.0, 1.0 
(0, 1.0 
(0, 1.0 
(0, 1.0 
(0, 0.75 
(0, 2 



Primary perturbation parameters in slow-roll parametrization 



Vera 

Vss 

( e ) 

ln(10"Ug) 

70 



Asz 



slow-roll parameter (when scale k = kg = 0.01 Mpc^^ exited the horizon); ryo-cr = 
slow-roll parameter (when scale fc = fco exited the horizon); r]a-s = g-L^ ^ , see 
slow-roll parameter (when scale k — ko exited the horizon); tjss = ^-L^ ^° ^ , see 
slow-roll parameter (when scale fc = fco exited the horizon); e — ig^c see 
Aq is the overall primordial perturbation power at fc = fco 



1 e„a„v 

8irG V 

Ref. |30] 
Ref. |30] 
Ref. 1301 



primordial ratio of correlated adiabatic component to total adiabatic power at fc 
[70 = sign(cos Ao) cos^ Ao, where Aq is the "correlation angle" at horizon exit of 
primordial isocurvature fraction at fc : 
amplitude of the SZ template 



= fco 
scale fco] 

fco; AqOo = primord. isocurvature power at fco 



(-0.075, 0.075 
(-0.075, 0.075 
(-0.075, 0.075 
(0, 0.075 
(1.0, 7.0 
(-1.0, 1.0 

(0, 1.0 

(0^ 



Ho 
h 

nar 
OcorO 



Derived parameters 
Hubble parameter [kms~^Mpc~^]; calculated from ojb, ojc, and 6 
h = iyo/(100kms-^Mpc-^) 
matter density parameter; Qm = (i^b + uJc)/h^ 
vacuum energy density parameter; flA = 1 — f2m 

spectral index of primordial uncorrelated adiabatic part; nar — 1 = dln(7'ar)/rflnfc|fe=fco 
spectral index of primordial correlated adiabatic part; riaa — 1 = dln(7'as)/dlnfc|fc=fe(, 
spectral index of primordial isocurvature part; niso — 1 = dln(Pg)/dln k\k=kn 
effective single adiabatic spectral index at fc = fco = 0.01 Mpc~^, Eq. (271 

Qcoro = sign(7o)^ao(l - ao)l7o|; ^o"corO = Cnsiko) = Cs7j(fco) = primord. correlation ampl. 



tophat (40, 100 
(0.40, 1.00 



total non-adiabatic contribution to the CMB temperature variance, Eq. ( 28 1 



Table I: Parameters and parametrizations. Our 11 (12 with primordial tensor perturbations) primary nested sampling 
parameters in two different parametrizations (amplitude and slow-roll) and a selection of derived parameters. In addition, in 
the amplitude parametrization, parameters Aq, 70, and ao can be derived from Ai, 7^, and Qi (i = 1,2). In the slow-roll 
parametrization, r is a derived parameter: ro = 16£(1 — j7oj). 



A. Phenomenological parametrizations 

The above-described model has six independent 
(scalar) perturbation parameters. We can choose the 
spectral index parametrization, where they are 

ItlAI, ao, 70, nar, »^as, ^iso, (29) 

or the amplitude parametrization, where they are 

InAf, \nAl ai, as, 71, \j2\ ■ (30) 



The model does not allow correlation to change sign as 
a function of scale, and therefore the sign of 72 is not an 
independent quantity. 

In addition, we can have primordial tensor perturba- 
tions whose power spectrum we write as 

Vrik) = Vriko) (fc/fco)"" • (31) 

We parametrize the primordial tensor power in the tra- 
ditional way by the tensor-to-scalar ratio 

Vrik) 



r{k) 



Vn{k) 



(32) 
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which may depend on the scale, as we see by substituting 
the power spectra into the above definition 



r{k) 



■To 



/ \ 



(l-|7ol) 1701(1;) 



where rg = r{ko). We assume the first consistency rela- 
tion (see, e.g., [HI [H [571 ESI [76]) 



'8(1 -|7|) 



(34) 



to be approximately valid at all scales of interest. This 
requires (moderate) "running" of tensor spectral index 
riT, since r and 7 depend on k. Hence we define ut up 
to first order as 



nrik) nriko) + ^qra \n{k/ko) , 
where the running is 

(f\nVT{k) _ dni 



qro — 



d(\nky 



din fc 



(35) 



(36) 



Demanding Eq. (34) to be valid at leading order in 
ln(fc/A:o) leads to 



qro = n-To [nra - {rim- - 1)] 



(37) 



This is the second consistency relation [76l I158j , and we 
impose this together with the first consistency relation as 
in '3] — see discussion after Eq. ( 42 ) on how Eqs. ( 34 ) and 



(371 follow from slow-roll inflation. Therefore, the pos- 



sible primordial tensor perturbations add only one extra 
parameter to the model, and we choose it to be tensor- 
to-scalar ratio ro at ko — 0.010 Mpc~^. This then gives 
the spectral index as a derived parameter, 



nro 



8(l-|7o|)' 



(38) 



and the running via Eq. (37). (In the phenomenological 



approach it would have been logical to keep r and as 
independent free parameters, and thus use the data to 
test the consistency relation. However, as neither ten- 
sor perturbations, nor deviations from adiabaticity have 
so far been detected, the data currently have not much 
power for this, and we choose to focus on the question of 
adiabaticity in this paper.) 

In phenomenological studies we choose the ampli- 
tude parametrization rather than the spectral index 
parametrization, since the amplitude parametrization 
leads to faster convergence of the parameter scan, leads 
to tighter and more reliable constraints on the parame- 
ters, and does not have the ambiguity of arbitrarily cho- 
sen prior ranges [T^ in Bayesian model selection that the 



spectral indices suffer from [1^1 [TT] . We list the primary 
(sampling) parameters of the amplitude parametrization 
in the first two blocks of Table [I] In the last block we 
list some interesting derived parameters, e.g., the spec- 
tral indices and amplitudes at fc = fco that are easy to 
calculate from the amplitudes at fci and fc2- 



B. Inflationary slow-roll parametrization 

In the slow-roll two-field inflation approach, we have 
the following relations between the slow-roll parameters 
and spectral indices and r, to first order in slow-roll pa- 
rameters [201 (see also, e.g., 1511 1571 1159] . and note 
that in Byrnes & Wands [30] n = stands for a scale- 
invariant spectrum whereas we have added the conven- 
tional 1 and we use a different sign convention for cos A 
and hence for tan A than [30] , since they have a different 
sign convention for TZ) 

"ar = 1 - 6e + 2?7o-a 

TT-as 1 - 2e + 2riss - 4rycrs tan A 

niso = 1 -2s + 2riss 

r = 16e sin^ A 

UT = -2e, (39) 

and to second order in slow-roll parameters, among oth- 
ers, 

qr - -8£2 + 4e?7„^ . (40) 
Here the correlation angle A is defined by 



cos A = — — — 7777 — sign[As-B)i 



(41) 



with < A < TT (so that sin A is non- negative) , and the 
relation to 7 is 

7 = sign(cos A) cos^ A 
sin^ A = 1 - I7I 
cos A = sign(7)y|7| 



tan A — sign(7)- 



(42) 



In principle, the slow-roll parameters are functions of 
scale fc, and the parameters we are fitting are the val- 
ues of the slow-roll parameters when our reference scale 
fco exited the horizon. We are using these values to con- 
struct the approximate forms of the primordial power 
spectra. 

In multi-field inflation, non-zero correlation (A ^ 7r/2, 
7^0) results in whenever the background trajectory in 
the field space is curved between the horizon exit (of the 
perturbations of cosmologically interesting scales) and 
the end of inflation. 
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From Eqs. ([39]) and Q we see that rg = 16e(l - |7o|), 
and riTo = — g(-xZ^'|.yp|) i ^-nd furthermore Eq. ( 40 ) can be 
rewritten as qro — — 2e[— 2£— (— 6£ + 2?7crcr)] = ?t-to['t-to~ 
('^ar ~ !)]• So, the tensor perturbations satisfy both the 
first and second consistency relations. 

After doing all our analysis using the phenomeno- 
logical amplitude parametrization we repeat everything 
using the theoretically motivated inflationary slow-roll 
parametrization. We summarize the parameters and no- 
tation in Table HI 



C. Specific examples 

Although we will compare to the data the two models 
described above, and uncorrelated or ±100% correlated 
variations of them, we review in this subsection three ex- 
amples on how a (correlated) mixture of adiabatic and 
isocurvature perturbations could arise without multiple 
(dynamical) inflaton fields. In these cases the spectral 
indices do not necessarily follow exactly the formulas of 
Eq. p9| . Neither are they fully free as in our phcnomcno- 
logical parametrization. So, constraining these specific 
models accurately would require additional parameter- 
scan runs. As a fourth example we review one of the 
simplest two-field inflation models, double quadratic in- 



where A = ("^x^) ' Since a is constrained to be 

small, the only observationally allowed possibility in this 
scenario is that A is small, which leads to (almost) un- 
correlated perturbations. 

If CDM is created from curvaton decay then the en- 
tropy perturbation is [160] 



S={l-R)2 



X* 



and 



aik) 



9(1 -Rf 
i?2(l + A-i) + 9(l-i?)2 



7(fc) 



1 + A 



(46) 



(47) 



flation, where the spectral indices obey Eq. (39) 



Now a can be small cither if A is small, which again leads 
to uncorrelated perturbations, or if R is very close to one, 
that is, if the curvaton dominates the energy density of 
the universe at the time of decay. In that case A can 
be arbitrarily large, leading to fully correlated pertur- 
bations, provided that the value of the curvaton field is 
sufficiently below Planck scale at the time when obser- 
vationally relevant scales exit the horizon. The spectral 
indices are given by 

War — 1 - 6£ + 27700 ; "as — "iso ~ 1 - 2£ , (48) 

where the parameters are evaluated at horizon exit. 



Mixed inflaton- curvaton scenario 



2. Modulated reheating with gravitino dark matter 



In many models the inflationary evolution is dominated 
by the inflaton field 0, but another light field is present 
during inflation and entropy perturbations in this field 
may be converted into curvature perturbations after in- 
flation is over. In the curvaton scenario the curvaton 
field X is constant during inflation but once the Hubble 
parameter becomes smaller than the curvaton mass, the 
curvaton begins oscillating and behaves like dust. Thus 
its energy density decreases slower than that of radiation 
and it can make a significant contribution to the energy 
density of the universe and to the curvature perturbation. 
The curvature perturbation is given by jl60| 



n 



1 



v/2Mj 



2R5x* 
3 X* 



(43) 



where R = 



is evaluated at the time of cur- 



.'^Px+'^Pi J 

vaton decay and "*" means that a quantity is evaluated 
at the time of horizon exit. 

If CDM is created before curvaton decay, the entropy 
perturbation is |160j 



S = -2R 



Sx* 



and we have 

a{k) = 



9A 

1 + lOA ' 



7(fc) = 



1 + A' 



(44) 



(45) 



In the modulated reheating scenario the decay rate 
r = r(cr) of the inflaton depends on an additional scalar 
field cr which was light during infiation. The inflaton de- 
cays into radiation when if ~ F, which means that due 
to entropy perturbations in the field a the inflaton will 
decay at different times in different parts of the universe. 
Assuming that after the end of inflation the inflaton oscil- 
lates and behaves like dust, its energy density decreases 
slower than that of radiation, which means that the parts 
of the universe where it decays later will have an over- 
density compared to the parts where it decays earlier. 
Thus the entropy perturbations in the modulating field 
are converted into curvature perturbations, given by 



1 



dT/da 



(49) 



If CDM is in the form of gravitinos then the entropy 
perturbation is given by 



1 fdV/da 

2 V F 



5(7.. 



(50) 



where + and — correspond to thermal and non-thermal 
production of gravitinos, respectively |161| . Then 



a{k) 



9g 
9 + lOC 



7(fc) = ± 



9 + e ' 



(51) 
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where 



(52) 



Again perturbations are constrained to be (almost) un- 
correlated. The spectral indices are the same as in the 
curvaton case. 



3. Axion 

In the axion model the axion field is massless during 
inflation, but obtains a mass after inflation ends, once the 
temperature of the universe drops below the QCD scale, 
Aqcd- When the Hubble parameter is smaller than the 
axion mass, the axion behaves like dust and can account 
for dark matter. The curvature perturbation in the axion 
model comes from the inflaton. 



n = 



1 



and the isocurvature perturbation from the axion, 

X* 



(53) 



(54) 



where w^/wc is the fraction of axions in CDM. The value 
of the axion field during inflation can be expressed as 
X* = fa^i, where fa is the axion decay constant and Oi 
is the initial misalignment angle. We now have 



= 7(fc) = 0, 



where 



A; 



Be* /w^ \ ^ / A/pi ^ ^ 



\^cj V fa 

The final abundance of axions is [H] 



-3 ^ f a'2 ^QCD \ I Ja 



(55) 



(56) 



= 10 X fd 



V200MeV 



IQiOGcV 



where /d < 1 is a dilution factor. The parameter fa is 
cosmologically constrained to be in the region 10^*^ GeV 
< /a < 4.1 xl0^2 GeV [13]. Thus if axions account for a 
significant portion of CDM, then e is constrained to be 
extremely small. The spectral indices are given by 



1 - 6£ + 2?700, 



1 - 2e. 



(57) 



4- Double quadratic inflation 



If there are several dynamically important fields during 
inflation then the direction of the background evolution 
plays the role of the inflaton. Perturbations in this direc- 
tion are adiabatic while perturbations perpendicular to 



it are entropy perturbations. The simplest such model is 
double quadratic inflation with two non-interacting fields 
with quadratic potentials. 



TA 1 2i.2|l 22 



(58) 



where > m^. Denoting the perturbations in the di- 
rection of the background trajectory by Sa and ones per- 
pendicular to it by Ss and assuming that the heavy field 
decays into CDM and the light field into radiation, the 
perturbations are |160j 



n 



S 



pi^ 



y2Af, 
4i?2 



sin 20*^* 



= (5(7. 



where tan 9 = x/4' siJ^d R 
y/(j)^ + R^X^. In this case 

a{k) = — 



(59) 

4e2(i?2_ 1)2 



where 



13 



/32 



32£, 
sin^ 29, 



7(fc) 



^-t-4e2(i?2 _ 1) 



4 f Mpi 



The isocurvature fraction a can be small either when 
/3 <C 1, in which case we can have arbitrary positive cor- 
relations, or when 2e{R^ — l) ^ (3 > 1, which means fully 
correlated perturbations^ The spectral indices for this 
model are given by Eq. (|39l with e = {d^Vf/ildnGV^), 



277s; 



7700(1 + cos 20) -t- 77p^x(l ~ cos 26*), 
7700(1 -cos26l)+77xx(l + cos 26*), (60) 



and r^crstanA — ijss, where ?7ij = {didjV)/ {SttGV). 



ABOUT THE ANALYSIS 



Data 



First we use only the CMB temperature and polariza- 
tion anisotropy data: Seven- Year Wilkinson Microwave 
Anisotropy Probe (WMAP7) data [H IT62] . the com- 
plete Arcminute Cosmology Bolometer Array Receiver 
(ACBAR) data J,, and QUEST at DASI (QUaD) data 
d] (QUEST = Q and U Extragalactic Survey Telescope, 
DASI = Degree Angular Scale Interferometer) . Then we 
add into the analysis either the supernova (SN) data from 
Supernova Cosmology Project (SCP) Union 2 compila- 
tion [4] or the large-scale structure data from Sloan Dig- 
ital Sky Survey Data Release 7 Luminous Red Galaxies 
(SDSS DR7 LRG) [5] that contains baryon acoustic os- 
cillations (BAG) [5]. 
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We perform the full Bayesian evidence calculation for 
each combination of data sets. Although this is compu- 
tationally demanding, it is important not to blindly com- 
bine different types of data sets without testing what is 
the individual information gained from each data set and 
whether they are consistent with each other. In particu- 
lar, the black-box method of combining all available data 
may lead to artificially tight constraints, if there is 'ten- 
sion' between the data sets, and there is between CMB 
and large-scale structure data. One historical example is 
using big-bang nucleosynthesis and observed abundances 
of light elements to constrain wt; tension between differ- 
ent elements led to too tight constraints, and a wrong 
value for (see discussion in, e.g., 1163,). Therefore we 
strongly advocate the approach of adding data sets to 
the analysis one by one. 

CMB: Our choice of the CMB data sets that extend 
the coverage from WMAP7 to smaller scales and also 
provide extra polarization data are ACBAR and QUaD. 
We impose in the CMB likelihood calculations a cut-off 
multipole -^max = 2100. We do not use at the South Pole 
Telescope (SPT) data 164 nor The Atacama Cosmology 
Telescope (ACT) data 165J that would extend with high 
precision to even much higher ^, for the following reasons: 
1) These data were not publicly available when we started 
our massive MultiNest runs. 2) In order to employ these 
data correctly the foreground modelling is of vital impor- 
tance, requiring marginalization over many "extra" nui- 
sance parameters that would slow down the convergence 
of our parameter scan. 3) Since the tensor spectrum in 
our analysis is red [ut < 0, see Eq. (34)], the tensor con- 
tribution modifies only the low-^ part of spectrum. So, 
the SPT and ACT data would not give extra constraints 
on tensor-to-scalar ratio. 4) As we will see, in the in- 
flationary parametrization also the isocurvaturc modifies 
only the low-i spectrum. The only effect of the SPT 
and ACT data would be possibly a slightly tighter con- 
straint on n-iso in the phenomenological parametrization, 
and somewhat tighter constraints on the background pa- 
rameters that could indirectly constrain nonadiabaticity. 

SN: There are many SN data sets on the market. 
Based on [166], our aim was to choose a recent set that 
would be the most "conservative" , giving the loosest up- 
per bound for the nonadiabatic contribution in the phe- 
nomenological mixed model that prefers larger than 
the adiabatic model [HI HZl US] . Since the SN data con- 
strain r^A, the largest isocurvature fraction is allowed 
by the data that allow for the largest ft a. Even the 
same "raw" SN data give very different results depend- 
ing on how the light curves are fitted. In general, the 
Spectral Adaptive Lightcurve Template (SALT-II) fitter 
[167) leads to a Hubble diagram that prefers a larger J7a 
than the Muhicolor Light Curve Shapes (MLCS2k2) fit- 
ter |168j in the ACDM model [166] . Therefore, our choice 
is SCP Union 2 compilation that has been obtained using 
SALT-II. In order to be as conservative as possible, we 
run the likelihood code that comes with CosmoMC with 
systematic errors turned on. 



MPK: It should be noticed that in Sec. lVIII we use the 
full matter power spectrum data (MPK) as a function of 
k that include the baryon wiggles spectrum BAO(A:), not 
the simplified BAO distance measure(s) or BAO scale(s). 
In the mixed models, these simplified measures cannot be 
used as given in literature, since the extraction of these 
numbers from the data assumes adiabatic initial condi- 
tions, see [24j (and \169 - 171 ) for detailed accounts on 
the danger of misinterpretation of these measures in the 
presence of isocurvature modes. Indeed, employing the 
MPK data with BAO in a right and consistent way in 
the mixed models is not a trivial task. The detailed 
procedure and modifications we did to CosmoMC and 
the MPK likelihood code [5] will be described elsewhere 
[172] . In short: for each parameter combination, we ex- 
tracted the oscillations separately for the adiabatic and 
isocurvature modes with unit primordial amplitude and 
scale-invariant spectrum by going from Fourier space to 
real space, i.e., by finding the correlation function. This 
has a bump that was removed with a technique described 
in the Appendix A.l of |173j . The bumpless correla- 
tion functions were transformed back to /c-space. These 
smooth adiabatic and isocurvature matter power spec- 
tra were then weighted by the initial amplitudes and the 
spectral tilts were added. Finally, the two smooth adi- 
abatic spectra Pai{k) and Pas(^) as well as the smooth 
Piso(fc) and sign(7)V^as^'iso (the correlation) were added 
together. The same was done for the original component 
spectra giving the unsmooth spectrum with the bary- 
onic wiggles, Pfuii(fc). After this the MPK likelihood 
routine was used in the normal way [S], i.e., a nonlin- 
earity correction was applied to the smooth spectrum, 
and, for presentation purposes, BAO was calculated as 
BAO(fc) = logio[PMi(fc)/Psmooth(fc)]- When using the 
MPK data, in order to reliably find Psmooth we restrict 
the prior ranges of ujb to (0.016, 0.029), uj^ to (0.05, 0.17), 
1006* to (0.90, 1.18), and r to (0.02, 0.20), while keeping 
the priors of all the other parameters the same as indi- 
cated in Table m 



B. Sampling method 

For the full parameter scan and for calculating the to- 
tal model probability Z (or Bayesian evidence) of each 
variant of our model as well as the reference adiabatic 
models we employ MultiNest package |26j that is based 
on an efficient variant of nested sampling algorithm. 

After the conceptual introduction by |174j . this 
method has been first applied to cosmology in a sim- 
ple case by |175j . Its most sensitive part, the sam- 
pling technique, has been refined by |176j and |177j 
to minimize the required number of likelihood eval- 
uations and to deal efficiently with possible patholo- 
gies, such as multi-modal posterior distributions and 
strongly curved parameter degeneracies. Finally, an even 
more robust and efficient code has been released by [22] 
for applications in cosmology, astronomy and particle 
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physics. The package, available for public use from 
|http://ccpforge .cse.rl.ac.uk /gf /pro j ect /mult inest / , con- 
tains an easily usable interface for the CAMB/CosmoMC 
cosmology code [T75J I179| that we have modified to han- 
dle an arbitrarily correlated mixture of adiabatic and 
isocurvature perturbations. 

The MultiNest user has simply to tune three parame- 
ters: the tolerance (accuracy) , the number of live points 
N, and the maximum efhciency e, which sets how aggres- 
sively we want the reduce the volume of parameter space 
at each iteration. A very attractive feature of MultiNest 
is that the need of a proposal matrix for the parameters' 
covariance, a well known hassle for conventional Markov 
Chain Monte Carlo (MCMC) users, is now superseded. 

In this paper we set the efficiency parameter in Multi- 
Nest to 0.3, the tolerance parameter to 0.5, and the num- 
ber of live points to iV = 400. With these settings the 
error estimate for the logarithm of the total model prob- 
ability turns out to be |(51nZ| ~ 0.20 - 0.26 for all of our 
runs. Therefore, when comparing model q to model q' , 
the logarithm of the Bayes factor, i.e., the logarithm of 
the ratio of the model probabilities, 

Bf = InBgy, = AlnZ = ln(Z,/Z,0 (61) 

has an uncertainty SBf ~ ±0.4. 



C. Bayesian model comparison and Jeffreys' scale 

A brief summary of Bayesian model comparison, cal- 
culation of evidence, i.e., the model probability P — Z, 
the sampling method, and MultiNest can be found in the 
Appendix of |17j . For a comprehensive review of many 
aspects of Bayesian cosmology, see |180j . 

The Bayes factors are classified with Jeffreys' scale 
[181) for the strength of evidence. This empirically cali- 
brated scale has thresholds at certain values of odds. If 



we take the model q in Eq. ( 61 ) as our null hypothesis and 



the model q' as an alternative model, then the original 
JeflFreys' book [ISlj grades verbally the odds according to 
Table im Jeffreys calls the simpler model (which is nested, 
i.e., a subset of the alternative model — with the extra 
parameters of the alternative model fixed, often to zero 
values) the "null hypothesis" , and the more complicated 
model that would lead to a new discovery (and deter- 
mination of extra parameters) an "alternative" model. 
From this set-up it is clear that one can only find evi- 
dence against the null hypothesis, not evidence against 
the alternative model, but merely just an indication that 
the null hypothesis is supported (while the extra param- 
eters of the alternative model could be so small as not to 
show up) in light of the current data. 

In the recent literature, q and g' have often been con- 
sidered from symmetrical footing, and indeed this is the 
only approach one can take in situations where q is not 
a subset of q' , i.e., whenever the two models are totally 
different and don't share (any) parameters. However, in 
all our cases we have a clear null hypothesis and hence 







Bf 


Interpretation 


> 1 




> 


"Null hypothesis q supported" 


1:3.16 


.. 1:1 


-1.15 


"Evidence against q, but not 








worth more than a bare mention" 


1:10 ... 


1:3.16 


-2.72 


"Evidence against q substantial" 


1:31.6 


.. 1:10 


-3.45 


"Evidence against q strong" 


1:100 . 


. 1:31 


-4.61 


"Evidence against q very strong" 


< 1 : 100 


< -4.61 


"Evidence against q decisive" 



Table II: A verbal interpretation of various odds (Bggi) 
and the corresponding logarithm of the Bayes factor {Bf = 
In Bgg' ) at the lower bound according to Jeffreys [181| . 





= Pq/Pq' 


Bf= In Bgg, 


Interpretation 


> 3 : 


1 


> +1.0 


"'Null hypothesis q supported" 


1:3 .. 


3:1 


-1.0 ... +1.0 


"Inconclusive" 


1:12 


.. 1:3 


-2.5 ... -1.0 


"Weak evidence against q" 


1:150 


... 1:12 


-5.0 ... -2.5 


"Moderate evidence against q" 


< 1 : 


150 


< -5.0 


"Strong evidence against g" 



Table III: Simplified thresholds and modernized verbal inter- 
pretation of various odds (Bggi) and the corresponding loga- 
rithm of the Bayes factor (Bf). The tabulated values of Bggi 
are rounded, while the values of Bf are the actual thresholds 
used in most of the modern cosmology literature. The actual 
wordings vary from source to source. 



we take a conservative "Jeffreys' approach" : When com- 
paring adiabatic and mixed models the null hypothesis 
is the adiabatic model, where the "nonadiabaticity" pa- 
rameters have null values. When comparing models with 
and without primordial tensor contribution, the null hy- 
pothesis is the model without tensors where Tq or e has a 
null value. In models with and without CMB lensing the 
number of parameters is the same, but the model with 
lensing has some extra ingredients, e.g., lensing potential 
that modifies the unlensed Ci spectrum, and in this sense 
it is natural to choose the model without lensing as our 
null hypothesis. For classifying the evidence against null 
hypothesis, we will use simplified thresholds and modern- 

in order to take 



ized wording |180j . modified in Table III 



into account our conservative null hypothesis approach 



D. Lensing: Bayesian evidence and a degeneracy 
with isocurvature 



All runs were made with CMB lensing on, unless stated 
otherwise. Initially, we did most of our runs without 
lensing, which is by a factor of 4-10 faster. Compar- 
ison of lensed and unlensed runs showed that lensing 
had only a minor effect on most of the 1-d posterior 
probability densities. However, the best-fit improve 
a lot when allowing for lensing (Ax^ ~ —8 in the adi- 
abatic case, —6 in the mixed case), and furthermore 
we find moderate Bayesian evidence against the unlensed 
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models. For example, in the amplitude parametrization 
the adiabatic unlensed model (without tensor perturba- 
tions) has InZ = —3861.3 while the lensed model has 
\nZ ~ —3856.6. This gives the for the logarithm of Bayes 
factor Bf = In Zunicnscd - In Zionscd = -4.7 ±0.4, or odds 
-Punicnscd/fionsod = exp(-4.7) = 1 : 110. In the mixed 
model the odds against the unlensed model are less over- 
whelming: Bf = -3.3 ±0.4, i.e., Punlcnsod/flcnsod = 1 : 

27. 

The reason for this is that in the amplitude 
parametrization large isocurvature spectral indices riiso ^ 
2 are favored and this causes the isocurvature and corre- 
lation contribution to affect rather uniformly the whole 
spectrum [21]. In particular, at high multipoles the al- 
most opposite phase isocurvature (and correlation) com- 
ponents are able to "fill" dips and "eat" peaks of the 
adiabatic spectrum, and hence they smear the total Ce 
spectrum compared to the pure adiabatic case. The CMB 
lensing has a similar smearing effect. Since the peaks 
and dips at high i in the data are smoother than the 
unlensed pure adiabatic ACDM model predicts, the fit 
to the CMB data improves considerably either by al- 
lowing for the isocurvature contribution (in the unlensed 
case Ax^ « — 5.5 for the best-fit mixed model in ampli- 
tude parametrization compared to the best-fit adiabatic 
model) or by turning the lensing on. Therefore, there is 
a degeneracy between lensing and isocurvature, and the 
Bayesian evidence against "no lensing" is weaker in the 
mixed model. 



VI. THE GENERAL CASE WITH CMB DATA 

A. Phenomenological approach — Amplitude 
parametrization 

An important question in constraining cosmologies is 
how much the assumptions made in the analysis af- 
fect the values of cosmological parameters obtained from 
the given data |182j . Therefore, we start by compar- 
ing the marginalized 1-d posterior probabilities in the 
mixed adiabatic and isocurvature model to the results 
obtained for the pure adiabatic model, using the ampli- 
tude parametrization. Fig. [l] shows these for the primary 
parameters and Fig. [2] for selected derived parameters 
for models with and without primordial tensor perturba- 
tions. In Appendix [a] we tabulate 68% or 95% confidence 
level (C.L.) intervals for selected parameters as well as 
Bayesian evidences Z for each model. 

We notice that the tensor contribution does not sig- 
nificantly affect any of these results, and the maximum 
allowed primordial tensor-to-scalar ratio, rp, is insensitive 
to the inclusion of the isocurvature mode: tq < 0.249 in 
the mixed model and ro < 0.264 in the pure adiabatic 
model at 95% C.L. The reason for the lack of interplay 
between tensor and isocurvature constraints is that the 
consistency relation requires the tensor perturbations to 
have a red spectrum (ht < 0), whereas the isocurvature 



perturbations favored in our phenomenological approach 
have a steeply blue spectrum. The slight formal tight- 
ening of the constraint in generally correlated models is 
due to the use of the first inflationary consistency relation 
that with a fixed ro leads to more negative ht (whenever 
7 7^ 0) than in the adiabatic model, see Eq. ([34]), and 
hence a larger effect on large scales. This has been no- 
ticed also in [H] . 

There is no evidence for the presence of tensor pertur- 
bations. Indeed, Bayesian model selection supports the 
model without a tensor contribution: in the mixed model 
the odds are P„o tensors /-Ptonsors = 7 : 1 (S/ = +1.96), 
and in the adiabatic model 11 : 3 « 3.7 : 1 {Bf = +1.30) 
in support of the model without tensors. Nevertheless 
we choose the model with primordial tensor perturba- 
tions as our baseline case for three reasons: 1) In infla- 
tionary models the tensor perturbations arise naturally, 
and hence it is more straightforward to compare our phe- 
nomenological model to the inflationary slow-roll model 
if we include the tensor perturbations. 2) Given point 
1, we have a theoretical prejudice that the tensor per- 
turbations exist, but probably with much smaller ampli- 
tude than we guessed when setting our linear flat prior 
ro < 0.75, which is an integral part of the "model" in the 
Bayesian evidence calculation and which we had to de- 
cide "before seeing the data" . So, what we found above 
is that the model without tensors (ro — 0) is supported 
over the model with ^ 0(0.3). 3) In the literature and 
in our previous work the model without tensor perturba- 
tions has already been studied extensively, though with 
older data sets than here. 

The CMB data allow a CDM (baryon) isocurvature 
fraction a < 0.064 (0.42) at scale ki, a < 0.154 (0.67) at 
scale fcg, and a < 0.512 (0.95) at scale ^2, but does not 
show any clear preference for a nonzero isocurvature con- 
tribution. The transfer function that maps the primor- 
dial entropy perturbation S into today's angular power 
spectrum Cf° has an extra factor of (keq/k)'^ oc {£eq/()'^ 
compared to the transfer function that maps the primor- 
dial curvature perturbation into C|'' and Cf^ |25l 11831 - 
I185j . (Here k^g is the wave number of perturbation that 
re-enters the Hubble horizon at matter-radiation equal- 
ity.) Therefore primordial isocurvature perturbations 
contribute much more strongly to large-scale (small fc, 
small £) CMB than to small-scale CMB, when compared 
to adiabatic perturbations of the same primordial mag- 
nitude, and hence the limits to a are tighter at large 
scales than at small scales. Indeed, if no (extra) prior 
is placed on the isocurvature spectral index, our phe- 
nomenological approach formally favors a large niso ~ 2 - 
3, i.e., a large primordial isocurvature fraction at small 
scales. For such isocurvature perturbations, a positive 
correlation with curvature perturbations is favored. The 
marginalized posterior probability density for the nona- 
diabatic contribution to the CMB temperature variance, 
ax, appears to show a weak detection of a positive value 
(ttT > at 90% C.L.), but Bayesian model compari- 
son doesn't support the inclusion of primordial CDM (or 
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Figure 1: General case, amplitude parametrization, 

CMB data. Marginalized 1-d posterior probability densities 
of the primary parameters (except Asz which would be flat) 
in amplitude parametrization. Red lines are for the purely 
adiabatic model, and black lines for the mixed adiabatic and 
isocurvature model. The solid lines are for the models with 
primordial tensor perturbations, while the dashed lines are 
for models with only scalar perturbations. 
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Figure 2: General case, amplitude parametrization, 

CMB data. Marginalized 1-d posterior probability densities 
of selected derived parameters from runs made in amplitude 
parametrization. The line styles are the same as in Fig. [l] 



baryon) isocurvature mode: PadiabaticZ-Pmixcd = 240 : 1 
{Bf = +5.48) with tensors, or 137 : 1 {Bf = +4.92) 
without tensors. The last column of the first two rows 
of Table |V] indicate that the CMB sets an upper limit 
of 5% for ax which is our pivot-scale free measure of 
nonadiabaticity. 
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Figure 3: General case, amplitude parametrization, 

CMB data. Marginalized 2-d posteriors of selected parame- 
ters for the mixed model with tensors. The thick black curves 
indicate the 68% and 95% C.L. regions. The intersection of 
thin black (red) lines is the best-fit mixed (adiabatic) model. 
The difference between these two best-fit models is Ax^ ~ —3 
in favour of the mixed model. The colour scale indicates the 
mean likelihood. 



According to Fig.[T]the inclusion of isocurvature affects 
significantly the sound horizon angle 6: mixed models 
prefer larger 9 than adiabatic models. They also prefer 
slightly smaller ujc- These changes are reflected in the 
derived parameters Hq and Qa in Fig. [2] mixed models 
prefer larger values. Fig. [3] demonstrates this degeneracy 
between the nonadiabatic contribution and 9 or 

The reason for the above degeneracy is extensively ex- 
plained in [9] and |17j . see in particular figure 2 and its 
explanation in |17) . In short: The amplitude and phase 
of the correlation component is somewhere in between 
the adiabatic and isocurvature components. Therefore 
the correlation component gives the main "nonadiabatic" 
contribution to . When the isocurvature and hence the 
correlation spectral indices are much larger than 1, the 
positive correlation component adds some power to the 
right end of each acoustic peak, moving them slightly 
to right and making them wider; see figure 2 in [3]. A 
larger 9 counteracts this effect by squeezing and moving 
the whole acoustic peak structure to left, hence retaining 
the almost "adiabatic" shape of the temperature angu- 
lar power spectrum that fits well the data. The sound 
horizon angle is defined hy 9 = Ts/Da, where Vg is the 
sound horizon at last scattering and Da is the angular 
diameter distance to last scattering. As depends only 
on LOh and lOc that are determined by the relative heights 
of the acoustic peaks, and as even this dependence is very 
mild |r86j . the CMB constraint on 9 is directly reflected 
in the favored DA^ which depends on f2,„, (or D,x in 
spatially curved models), and Hq. In the curved mixed 
models a well-fitting Da is achieved in a slightly closed 
case [T7], i.e., when il,„ + ilA ^ 1.05, but in the flat mixed 
models studied here, the constraint + ~ 1 dictates 
that the only way to achieve a well-fitting value for Da 
is to increase JIa (and Hq) and decrease correspondingly 
il,„ compared to the well-fitting adiabatic case. 

Forcing the background parameters, for example f^A, 
closer to the "adiabatic" values should set indirect extra 
constraints on isocurvature or at least disfavour the pos- 
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itive correlation and large niso- Since fl^ is constrained 
by the supernova data and ftm (= 1 — ^^a) by the large- 
scale structure (matter power spectrum) data, we check 
in Sec. |VII| how our results are affected by the inclusion 
of these data. However, it should be emphasized that the 
extra constraints on J^a may lead to tighter constraints 
on isocurvature only in the spatially flat case which we 
consider here, leaving the constraints untouched if spa- 
tially curved geometry is allowed [T7]- 

In addition, the amplitude and phase of baryonic wig- 
gles (BAO) in the matter power spectrum might set di- 
rect extra constraints on the isocurvature and correla- 
tion components if their relative amplitudes at scales 
probed by MPK were large enough. As the CMB tightly 
constrains the isocurvature contribution on large and 
medium scales (wave numbers fci and fco, respectively), 
the potential extra constraints from BAO come in cases 
with a large nonadiabatic component on small scales 
(k ^ ^2), i-e., whenever riiso, and hence n^or, are very 
large, riiso 2.5. However, we will show in Sec. VH that 



we would need more precise data than the current SDSS 
DR7 LRG data. 



Two-field inflation approach 
parametrization 



Slow-roll 



Fig. [4] shows marginalized 1-d posterior probability 
densities for the primary parameters in the slow-roll 
parametrization, and Fig J5] selected derived parameters. 
We have collected in Fig. [6|some 2-d combinations of the 
four slow-roll parameters in the mixed model with ten- 
sors. As indicated in Table |l| our prior constraint is that 
the magnitudes of the slow-roll parameters are smaller 
than 0.075, i.e., <C 1. We see that the data are not able 
to constrain the slow-roll parameters this well: the 95% 
C.L. contours intersect this prior constraint, except for e, 
which we constrain to £ < 0.048 at 95% C.L. This is three 
times weaker than the limit e < 0.016 in the adiabatic 
model, see Table |V] in Appendix |A] Using wider priors 
would have the problem, that our approximation of cal- 
culating the spectra only to first order in slow-roll param- 
eters (see Sec. IV B I would become questionable; includ- 



ing higher-order terms, on the other hand, would require 
inclusion of additional slow-roll parameters, whose de- 
termination would be infeasible. Thus the current CMB 
data are not yet sufficiently accurate for usefully con- 
straining the slow-roll parameters of two-field infiation 
models, except for s. 

The assumption that the magnitudes of slow-roll pa- 
rameters are small, requires that all spectra are close to 
scale invariant [see Eq. ( 39 )] , except the correlated part of 
the primordial curvature perturbation in the cases where 
7o is close to zero.^ These cases are excluded by the 



We see from Eq. (|23| that the relative amplitudes of the un- 



data, so in practice also rXas will be close to 1. In par- 
ticular, the isocurvature perturbations are not allowed 
to have a very different spectral index from the adiabatic 
perturbations. This means, that the isocurvature pertur- 
bations contribute much more to large-scale CMB than 
to small-scale CMB, and the constraints to the isocurva- 
ture fraction come mainly from the comparison between 
the level of large- and small-scale CMB anisotropy, and 
not from the location of the acoustic peaks in the angular 
power spectrum. The near scale invariance of all spectra 
also means that the upper limit to a is roughly the same 
at all scales. The upper limit ao < 0.026 turns out to be 
even tighter than our large-scale {k = 0.002 Mpc~^) up- 
per limit ai < 0.064 in the phenomenological approach. 

Since tensor perturbations share with isocurvature per- 
turbations the property that, when close to scale invari- 
ant, they contribute more to large-scale CMB than to 
small-scale CMB, we now get a degeneracy between the 
tensor and nonadiabatic contribution. In Fig. [7] we show 
how the allowed tensor contribution, parametrized by 
the primary parameter e or the derived parameter ro, 
depends on the nonadiabaticity parameters 70, ao, and 
Q!corO- Because negative correlations give a negative con- 
tribution to large-scale C^, whereas tensor perturbations 
give a positive contribution, a larger tensor contribution 
is allowed when the isocurvature mode is negatively cor- 
related than when it is positively correlated, see the left 
hand side panels of Fig. [7] Or the other way around, 
when allowing for the primordial tensor contribution, 
the 95% interval of Ofcoro changes from the symmetri- 
cal (-0.092, 0.092) to (-0.101, 0.092) while the median 
moves from -0.021 to -0.034, see Table [V] However, since 
strongly correlated isocurvature perturbations cause a 
much larger signal for a given ap than uncorrelated ones, 
the limit to the CDM isocurvature fraction is formally 
tighter (ao < 0.026), when tensor perturbations are al- 
lowed, than when they are not (ao < 0.068). The cor- 
responding numbers for the baryon isocurvature are 0.29 
and 0.44. 

Allowing for isocurvature perturbations also leads for- 
mally to a tighter limit on the amplitude of tensor pertur- 
bations (ro < 0.18) than in the adiabatic case (ro < 0.26). 
The main reason for this is the tan A term in the expres- 



correlated and correlated adiabatic components are 1 — I70I and 
I70I at scale fco- Assuming the slow-roll parameters to be close 
to zero gives, according to Eq. (jSOjl, riar = 1 — 6e + 2riacr 1 
and rias = I — 2e + 2r]ss — '^Virs tan A ~ 1 — Arja-s tan A. Now 
the relative amplitudes at another scale, fc, are 1 — I70I and 
|7o|(fc/fco)"''''''=*'""^- Assuming 70 = -0.0001 and jj^s = 0.01 
we find tan A —100 and —4r]as tan A 4. In this case, at 
k > 8.4fco (or at fc < 0.12A;o if r)<js x 7o had the opposite sign), 
the correlated adiabatic component would be the dominant com- 
ponent even though, at first sight, we would naively think that 
I70I = 0.0001 means a negligible contribution from "Pas- This ex- 
ample shows that there are cases where an extremely small I70I 
leads to a large rtas and to a large Pas either at the large-scale 
end or the small-scale end of the CMB angular power spectrum. 
However, these models are ruled out by the data. 
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Figure 4: General case, slow-roll parametrization, 

CMB data. Marginalized 1-d posterior probability den- 
sities of the primary parameters (except Asz) in slow-roll 
parametrization. The line styles are the same as in Fig. [l] 
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Figure 5: General case, slow-roll parametrization, 

CMB data. Marginalized 1-d posterior probability densities 
of selected derived parameters from runs made in slow-roll 
parametrization. The line styles are the same as in Fig. [l] 



sion for rias in Eq. ( 39 ) . As explained in footnote [T] if 70 
was too close to zero, the V^s spectrum could be steeply 
red or blue, and hence excluded by the data. This disfa- 
vors those values of 70 that would lead to a large tq, since 



according to Eqs. ([39) and Q = 16£(1 - I70I). More- 
over, due to the 1^7o| factor, and due to large (nega- 
tive) correlations being favored, the constraint e < 0.048 
is much weaker than in the adiabatic case (e < 0.016). 

The 1 — I70I factor in tq together with the upper left 
corner plot of Fig. [7] also explains the difference of 1-d 
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Figure 6: General case, slow-roll parametrization, 
CMB data. Marginalized 2-d posteriors for the two-field 
inflation slow-roll parameters from MultiNest runs with pri- 
mordial tensor contribution allowed (e > 0). 
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Figure 7: General case, slow-roll parametrization, 

CMB data. Marginalized 2-d posteriors for various com- 
binations of nonadiabaticity (70, ao, acoro) and tensor per- 
turbation parameters (e, vq). 



posteriors of 70 in Fig. [4j The data constrain the pri- 
mordial tensor amplitude, tq, but whenever I70I is close 
to 1, we obtain an acceptably small value for tq no mat- 
ter what value the primary parameter e has. Therefore, 
when marginalizing over s, there is a large volume of 
well-fitting models when I70I is close to 1. In contrast, 
when 7o is close to zero, e must be small in order to ob- 
tain a small Tq. Hence the well- fitting parameter space 
volume is small, if I70I is small. This means that the dif- 
ference between 1-d posteriors of 70 of "no tensors" and 
"with tensors" cases is due to a marginalization artifact, 
and the apparent preference of 70 = —1 with tensors 
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is not strictly speaking physical. From the mean like- 
hhoods (indicated by the colour scale) in Fig. [t] we see 
that indeed the best-fitting models have only a mildly 
negative correlation, 70 ~ —0.1, and the good- fit region 
covers the range (—0.95, 0), but does not particularly 
favour 7o — —1. Indeed, often 1-d marginalized poste- 
riors give limited information on high-dimensional prob- 
lems, in particular, if some parameters are loosely con- 
strained: although the 1-d Figs. |4] and [5] give only weak 
hints, the 2-d Fig. [7] supports the conclusion that in the 
slow-roll parametrization, negative correlation (negative 
7o and acoro) is clearly favoured — slightly more favoured 
with tensors than without them. 

The amplitude parametrization favours a positive cor- 
relation. So why does the slow-roll parametrization 
favour a negative correlation? The explanation lies again 
in the spectral indices, i.e., what part of the angular 
power spectrum the nonadiabatic contribution is able to 
modify. In the amplitude parametrization the width (and 
location) of the acoustic peaks got slightly modified due 
to large njgo and ricor- In the slow- roll parametrization 
the isocurvature mode does not have the ability to mod- 
ify this region, since riiso ^ n-ar- So, the nonadiabatic 
contribution can significantly modify only the large-scale 
{i < 50) part of C^. There the angular power at some 
multipoles is lower in the data than predicted by the 
adiabatic ACDM model, most notably the quadrupole, 
£ = 2. An uncorrelated or positively correlated isocurva- 
ture contribution is thus disfavoured as they would fur- 
ther increase the discrepancy between the data and the- 
ory, but a negative correlation can reduce the power in 
this region. However, due to the cosmic variance the con- 
tribution to from this region does not have as much 
weight as the one from the acoustic peak region, and thus 
the difference of the best-fit mixed model and adiabatic 
model is only Ax^ w — 2 in favour of the mixed model 
(while it is —3 in the amplitude parametrization). 

The shapes of the preferred regions in the slow-roll 
parameter space can be understood by the preference for 
certain values of the spectral indices. From Eq. (39) we 
have that riar 1 corresponds to £ ~ rja-a/S. Since n^r 
somewhat less than 1 is preferred by the data, this line 
is shifted a bit to the left in the {rj^^, e) plane, see Fig. [6] 
The slope of the good-fit region in the (rissiVtrs) plane is 
about —2/3, i.e., 2riss — Srja-s ~ const., which corresponds 
to tan A w —3/4 or 70 ~ —0.64, which is in line with 
negative correlations being preferred. 

The effect of the optical depth r on the CMB temper- 
ature angular power spectrum is also somewhat degener- 
ate with isocurvature and tensor perturbations, affecting 
the relative level of low and high £, but since t is mainly 
determined by low-£ WMAP7 polarization data, where 
it causes the distinctive reionization bump, isocurvature 
perturbations have negligible effect on the determination 
of r. 

Finally, the odds for the adiabatic model compared 
to the mixed model in the slow-roll parametrization 
are 24 : 1 (6/ = +3.19) with tensors, and 80 : 1 
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Figure 8: General case, comparison of amplitude and 
slow-roll parametrizations, CMB data. Marginalized 1- 
d posterior probability densities of selected parameters from 
runs made in amplitude and slow-roll parametrizations, al- 
lowing for a tensor contribution. The solid lines are for the 
amplitude parametrization and the dot-dashed for the slow- 
roll parametrization, red color for the adiabatic and black for 
the mixed model. 



{Bf = -1-4.38) without tensors, in support of the adiabatic 
model. In the adiabatic model, with our priors for the 
slow-roll parameters, we find again support for the model 
without tensor contribution, Pno tensors /^tensors = 9:1 
{Bf = -1-2.23). In the mixed models this changes to 
Pno tcnsors/-Ptcnsors = 2.8 : 1 {Bf = +1.04), but recall- 
ing our error estimate SBf « ±0.4, the last result could 
fall into the category "inconclusive" . 



C. Comparison of parametrizations, and priors of 
derived parameters 

To allow for an easy comparison between our phe- 
nomenological and two-field inflation approaches, we re- 
peat the 1-d posteriors of some parameters of the adia- 
batic and mixed models with tensors in Fig. [Sj First we 
notice that the parameters of the adiabatic model are so 
well constrained that the two parametrizations lead to 
almost identical results. 

Since the allowed isocurvature contribution is smaller 
in the two-field inflation approach than in the phe- 
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Figure 9: Priors (dashed) and posteriors (solid) of selected 
derived parameters, and posteriors divided by priors (dot- 
dashed). The red color is for the amplitude parametrization 
and the black for the slow-roll parametrization. 



models with riiso ^ 3.5, see Fig. [9j and the slow-roll 
parametrization limits riiso close to one: 1 — 2 max(e) -|- 
2imii{riss) < niso < 1 — 2min(e) -t- 2max(775s). 

It is also interesting to compare how the parametriza- 
tion and prior ranges (that are considered as part of the 
model in Bayesian thinking) affect the Bayesian evidence. 
In the adiabatic case with tensors we find Bf — 1.30 and 
in the mixed case Bf = 3.59 in favour of the slow-roll 
parametrization. This is what one would have expected, 
since the prior ranges are so much narrower in the slow- 
roll parametrization. 

To complete the comparison of parametrizations, we 
show the temperature Ce spectra of two example models 
in Fig. |10[ Model I is from our amplitude parametriza- 
tion run and has roughly the maximum 95% C.L. allowed 
primordial nonadiabatic contribution at scale ^2 (acor2 = 
0.42). Model II is from our slow-roll parametrization run 
and has roughly the maximum allowed negative nona- 
diabatic contribution at ki (ofcori = —0.11). These 
models represent the above-described typical features 
of well-fitting mixed models in amplitude and slow-roll 
parametrizations, respectively. The parameters of these 
examples can be found in Tables [TV] and [V| in Appendix [5] 
under names "Model I" and "Model 11" . 



nomenological approach, the preferred values for the 
background parameters in the mixed model are closer to 
the adiabatic case in the slow-roll parametrization (see 
Fig. [8| . Even more important is the preferred negative 
sign of the correlation. Due to this, the preferred value of 
the sound horizon angle moves to the other side (smaller) 
of the adiabatic result, which is then reflected in Wc, Hq, 
and ^Ia- Since the negative correlation reduces the power 
at large scales compared to the small scales, the (slow- 
roll) models with negative correlation prefer a smaller 
spectral index n?^. 

It should be kept in mind that for the derived param- 
eters the prior is not necessarily flat, and hence the pos- 
teriors could in some cases directly reflect the prior or be 
biased by it. We have checked that this is not the case 
for our derived parameters. In Fig. [9] we show the prior 
in both parametrizations for selected derived parameters. 
We plot also the posterior divided by the prior, which in- 
dicates how the results would look like if we had a flat 
prior for the derived parameters. Importantly, the prior 
is almost flat for all of our derived parameters (including 
the isocurvature parameters) in the region of the peak of 
their posterior. There is one exception: since the CMB 
data do not constrain riiso, a flat wide range prior for it 
(like in the spectral index parametrization) would lead 
to very large values being "preferred" [El [25] , riiso ~ 3- 
7 depending on the chosen pivot scale. (See flgure 21 
in and flgure 5 in [T7]. Note from the latter fig- 
ure also that although ax is a pivot-scale free measure 
of non-adiabaticity, its posterior probability density is 
affected by the choice of pivot scale.). The amplitude 
parametrization leads naturally to a prior that excludes 



VII. THE GENERAL CASE WITH 
ADDITIONAL DATA 

We study how the results are affected by using super- 
nova (SN) [3] or matter power spectrum (MPK) [S] data 
in addition to the same CMB data as in the previous sec- 
tions. For brevity and clarity, we show figures only for 
models with tensors. Fig. [Tl] shows the 1-d posteriors of 
the primary parameters of amplitude parametrization for 
different data sets: 1) with CMB only, 2) with CMB and 
SN and 3) with CMB and MPK. Fig. [T2| gives the same 
sets for derived parameters. The analysis is repeated in 
the slow-roll parametrization in Figs. [13] and [T4j 

In the phenomenological approach, the additional data 
exclude those mixed models that have a large fl\ (see 
Fig. 12). These models have a large positive nonadi- 



abatic contribution, which thus gets less favoured (see 
7o, ttcorO: and ax)- Therefore, by constraining the back- 
ground parameters (see also Hq in Fig. 12 and Wc, in 
Fig. 11) closer to the "adiabatic values" or to less fa- 



vorable values to isocurvature, the positive nonadiabatic 
contribution gets indirectly constrained tighter than with 
the CMB data alone: the upper limit for ax changes from 
4.9% (CMB) to 4.2% (CMB and SN) or to 2.4% (CMB 
and MPK). 

The additional data have less effect in the slow-roll 
approach, since the background parameters were already 
closer to the adiabatic case. Actually the additional data 
prefer a smaller JIa than the CMB data in the adiabatic 
model, which is the same direction where the isocurva- 
ture contribution pulls in the slow-roll approach. There- 
fore, including the additional data make the mixed model 
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Figure 10: Typical temperature angular power spectra of well-fitting mixed models with a non-negligible nonadiabatic contri- 
bution. Model I (qt ~ -1-0.025) is from our amplitude parametrization runs and Model II (qt ~ —0.049) from our slow-roll 
parametrization runs. Both models are within Ax^ ~ 4 from the best-fit adiabatic model. Solid black curves are for the total 
lensed £{£ + l)Ce/{2n), the solid red for the adiabatic contribution, the solid green for the positive nonadiabatic contribution, 
the dashed green for the negative nonadiabatic contribution, and the dotted black for the tensor contribution. Parameters of 
these models are reported in Tables [IV] and [V] in Appendix [X] 



more favourable than with the CMB data alone. The up- 
per limit to the CDM isocurvature fraction relaxes from 
ao < 0.026 (CMB) to 0.032 (CMB and SN) or to 0.034 
(CMB and MPK). From the ar plot in Fig. [w] it is clear 
that the preference for small leads to a very strong 
preference for a negative nonadiabatic contribution. The 
95% C.L. range for ax in the slow-roll parametriza- 
tion changes from ar e (-5.8%, 4.5%) with CMB to 
(-5.6%, 3.0%) with CMB and SN or to (-6.3%, -0.8%) 
with CMB and MPK. So, in the slow-roll parametriza- 
tion the MPK data exclude positive nonadiabatic contri- 
bution with more than 95% C.L. 

The results with SN data are between the results of 
CMB alone and the ones with CMB and MPK data for 
most of the parameters. This happens, since the SN con- 
strain only mildly and are in a relatively good agree- 
ment with the CMB. In contrast, the MPK data prefer a 
significantly smaller (larger ftm) than the CMB OE]. 
If this discrepancy in the determination of Qm persists 
in the future, then MPK data will continue to have a 
dramatic indirect effect on the determination of isocur- 
vature. 

In the amplitude parametrization the MPK data might 
lead to also a direct extra constraint on nonadiabatic per- 
turbations. At first sight, the niso plot of Fig. 12 seems to 
support this: the MPK data exclude the largest isocur- 
vature spectral indices which the CMB data are not able 
to constrain. More clearly this effect can be seen in the 
Q!2 plot of Fig. [11] However, as we will see below, in 
practice the extra constraint on a2 comes almost fully 
from the background constraint on $7^. This can be also 
realized by noticing that the SN data (that are purely 
background data) have a similar, but milder, effect. 

The MPK constraints on nonadiabaticity in the slow- 
roll parametrization are (and will always be) purely in- 
direct background constraints, since the nearly scale- 



invariant isocurvature perturbations with the tiny am- 
plitude allowed by the large-scale CMB are not able to 
modify the matter power spectrum at all. In contrast, in 
our phenomenological approach, which (after the CMB 
constraints) allows large primordial isocurvature and cor- 
relation components on small scales, the slope of the mat- 
ter power spectrum at large k can be affected [ini US] by 
the correlation component (and by the isocurvature com- 
ponent on very small scales) ^. However, here the MPK 
likelihood code tries to soften the effect of the overall 
shape (by marginalizing over two nuisance parameters) as 
the exact shape is not known due to the nonlinear effects. 
Therefore we would not expect tight extra constraints 
from the slightly modified slope, but instead from the 
phase and amplitude of the baryon acoustic oscillations 
(BAO). Unfortunately, primordial isocurvature perturba- 
tions of the same amplitude as the adiabatic ones lead to 
20 - 150 times smaller matter power spectrum today over 
the observable scales. This is clearly seen in Fig. [T5|^a) 
where we plot the resulting matter power spectra from 
scale-invariant primordial spectra with primordial ampli- 
tude 1. The background parameters in Fig. 15 a) are the 
same as in Model I. 



In Fig. 15 b) we show a break down of the total mat- 
ter power spectrum into its components for the mixed 
Model I of Tables IIVI and Ivl This model has a maxi- 



^ Our phenomenological model with a free nigo differs from, e.g., 
Ref. ilSTj where the adiabatic and isocurvature components 
share the same spectral index — the situation that is close to 
our two-field slow-roll inflation approach. As the CMB data 
prefer predominantly adiabatic nearly scale-invariant perturba- 
tions, the common spectral index is forced near to 1, which in 
|187l leads to a conclusion that the isocurvature would not affect 
the matter power spectrum. This is in agreement with what we 
concluded in the case of slow-roll parametrization. 
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Figure 11; General case, amplitude parametrization, 
comparison of CMB, CMB and SN, and CMB and 
MPK data. Marginalized 1-d posterior probability densities 
of the primary parameters in amplitude parametrization with 
different data sets: CMB (soUd hne style), CMB and SN (dot- 
dashed), or CMB and MPK (dashed). The red color is for the 
adiabatic and black for the mixed model. 
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Figure 12: General case, amplitude parametrization, 
comparison of CMB, CMB and SN, and CMB and 
MPK data. Marginalized 1-d posterior probability densities 
of selected derived parameters from runs made in amplitude 
parametrization with different data sets. The line styles are 
the same as in Fig. [TTj 



Figure 13: General case, slow-roll parametrization, 
comparison of CMB, CMB and SN, and CMB and 
MPK data. Marginalized 1-d posterior probability densities 
of the primary parameters in slow-roll parametrization with 
different data sets. The line styles are the same as in Fig. |11| 
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[■ ^^i' ^i*^ tensors, with CMB and SN data 
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^"^"^^;a-j:^7^^'7-~"i — • ^ Mixed, with tensors, (with CMB data) 
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Figure 14: General case, slow-roll parametrization, 
comparison of CMB, CMB and SN, and CMB and 
MPK data. Marginalized 1-d posterior probability densities 
of selected derived parameters from runs made in slow-roll 
parametrization with different data sets. The line styles are 
the same as in Fig. |11[ 



mum allowed (at 95% C.L. by CMB alone) primordial 
non-adiabatic contribution at scale k — k2 ^ 0.05 Mpc^^ 
that falls in the middle of the log(fc)-range of the MPK 
data: acor2 = 0.42, 72 = 0.73, 02 = 0.43, and niso = 3.3. 
We choose this model in order to study why the MPK 
data rule it out. For this model the isocurvature (correla- 



tion) contribution to the matter power at ^2 is 120 (6.7) 
times smaller than the sum of the two adiabatic com- 
ponents. As seen in Fig. [Tsj^b), the correlation compo- 
nent makes the matter power spectrum to decrease less 
steeply compared to the adiabatic spectrum at highest 
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Figure 15: (a) Linear matter power spectra at redshift z = Q (today) resulting from scale-invariant primordial spectra with 
unit amplitude. The dashed red curve is for the adiabatic primordial perturbations (Af = \, — 0, on = 0, riar ~ 1), the 
dotted blue for the isocurvature {A^ = 1, 7i = 0, = 1, ruso ~ 1), and the green for their correlation. The corr espon ding 
thin solid curves are the smooth spectra from where the baryonic wiggles have been removed as described in Sec. VA (b) 



Linear matter power spectra at redshift z — for Model I of Tables |IV| and |V] in Appendix |A] The gray dots with error bars 
are the SDSS DR7 LRG data j5j as supplied with CosmoMC, and the dashed thin black curve is the linear power spectrum 
('tot') modified to take into account the BAO damping, the uncertainty in the non-linear modelling, and finally convolved with 
the observational window function (see equation (5) in [5]) in order to be compared directly to the data. The solid thin black 
curve is the smooth part of this. Note that plot (d) provides a "zoom-up" of the difference between these dashed (wiggly) and 
solid (smooth) thin black curves that are not easily visible in this plot, (c) The linear baryon acoustic oscillations of plot (a); 
the full oscillating curves divided by the smooth thin solid curves, (d) The baryon acoustic oscillations of plot (b): the solid 
black, red and green curves are the linear BAO(fc) at redshift z = Q for the total, adiabatic and non-adiabatic components 
of Model I, respectively. The corresponding dashed curves that have smaller amplitude are obtained from the linear curves 
at three different redshifts modified to take into account the BAO damping, the uncertainty in the non-linear modelling, and 
finally convolved with the observational window function. For comparison, the dashed thick gray curve indicates BAO(fc) for 
the best-fit (to CMB and MPK) pure adiabatic model. The gray dots with error bars are the SDSS DR7 LRG data divided 
by the smooth part of the best-fit adiabatic model. The first two data points [see plot (b)] fall outside of the scale of this plot, 
the first being 0.22 and the second one 0.10. It should be noticed that plot (d) is only for illustration purposes, and the actual 
MPK likelihood is calculated comparing the theoretical convolved matter power to the data, as shown in plot (b), and taking 
into account the full covariance matrix of the MPK data that accounts for the strong correlations between the data points. 



wave numbers probed by the data. (Note that the obser- 
vational window functions related to each data point are 
rather broad, extending well beyond the last data point.) 



From Fig. 15 [c) we see that, although the amplitude 
of the isocurvature BAO is damped compared to the adi- 
abatic BAO, the isocurvature would distinctively modify 
the phase of BAO if the isocurvature contribution in the 
matter power spectrum was of the same order or larger 



than the adiabatic contribution. As seen above, after 
the CMB constraints, the maximum non-adiabatic con- 
tribution is ~ 10%. Moreover, this comes mostly from 
the correlation component whose phase is rather close to 
the adiabatic component. According to Fig. [Tsjjd), the 
maximum non-adiabatic contribution to the total BAO 
is < 5% (compare the green and black curves). This is 
not detectable by the SDSS DR7 data. To set direct ex- 
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tra constraints on isocurvature, the matter power data 
should be about an order of magnitude more accurate, 
which could be within the reach of Euclid |188) . Anyway 
the future LSS surveys will be important in constrain- 
ing the isocurvature modes by breaking degeneracies that 
the CMB leaves between the background parameters and 
nonadiabaticity. ^ 

As the non-adiabatic contribution to BAO (A;) does not 
explain the very bad MPK likelihood of Model I, the 
main reason has to be the wrong turn-over scale [com- 
pare the position of the peak of the best-fit (to CMB 
and MPK) adiabatic matter power spectrum to that one 
of Model I in Fig. llSFb)]. In the ACDM model the 



turn-over scale is determined by the scale that enters 
the Hubble horizon at the time when the matter and 
radiation energy densities are equal. This scale is /cgg = 
LOm X 0.0729 Mpc~\ i.e., Kq/h = hn„, x 0.0729 Mpc"^ 
For Model I, that has f^,„ = 0.203 and h = 0.783, this 
gives keq/h — 0.0116 Mpc~^ whereas for the best-fit adi- 
abatic model we have ^l^n = 0.293 and h = 0.688 leading 
to keq/h = 0.0147Mpc~^ in agreement with Fig. 15 ^b). 



The rest of the bad MPK likelihood of Model I comes 
from the wrong phase of (even the adiabatic component 
of) BAO caused by the "non-adiabatic" values of the 
background parameters, i.e., the very small flm- 

To test the above observations further we report sev- 
eral MPK x^: The best-fit adiabatic model to CMB and 
MPK has an MPK = 47.0 (60.0) with (without) bary- 
onic wigg les. Model I has = 60.5 (67.8) with (without) 
wiggles. Furthermore, Model I has = 60.8 without 
non-adiabatic wiggles. This is about the same as with 
all the wiggles. So, as expected, the "wrong" phase of 
(the less than 5% contribution of) non-adiabatic wigg- 
less does not affect x^. Interestingly, dropping the non- 
adiabatic contributions from Model I gives = 58.2 
(62.3) with (without) wiggless, i.e., improves the fit only 
by Ax^ — —2.3 (—5.5). This gives an indication of 
the magnitude of the effect from the modified slope by 
the non-adiabatic contribution. Finally, changing the 4 
backgound parameters to the values of the best-fit adi- 
abatic model, the mixed model gives = 48.8. This 
is only 1.8 worse than for the best-fit adiabatic model, 
and shows clearly that the bad MPK x^ of high isocur- 
vature fraction mixed models come almost solely from 
the unfavourable values of the background parameters 
while only Ax^ ~ 1-2 comes from the too gradual slope 
and in practise nothing from the phase information of the 
non-adiabatic BAO component. 



^ Detailed forecasts for combining the near future Planck '189','190' 
CMB data with the predicted LSS data from Baryon Oscilla- 
tion Spectroscopy Survey (BOSS, see e.g. |191l 11921 ') and Ad- 
vanced Dark Energy Physics Telescope (ADEPT, see e.g. 11931 ) 
or Euclid-like survey are provided in [24) and |169) . Using the 
LSS data together with Planck is expected to bring down the 
uncertainty of the isocurvature fraction from a few percentage 
points to better than 1%. 



VIII. SPECIAL CASES WITH CMB DATA 

A. Uncorrelated perturbations 

In the uncorrelated case there is no V^sik) spectrum, 
and therefore we have only 9 independent parameters: 
the usual 4 background parameters 



and 5 perturbation parameters 

In Al , riar , ao , ^iso , ro , 

where 

"ar = 1 - 6e + 2r/o.o- 

"iso = 1 - 2£ -f 2?7ss 

ro = 16£ . 



(62) 
(63) 



(64) 



In the amplitude parametrization, the primary pertur- 
bation parameters are 

In Al, \nAl, ai , a2 , ro , (65) 

and in the slow-roll parametrization 

In^Q , ao , 77o-(T , Vss £■ (66) 

The marginalized 1-d posteriors for this model with the 
CMB data are indicated by the solid blue 7 = curves 
in Fi g. |16| for the amplitude parametrization, and in 
Figs. |17| and [T8| for the slow-roll parametrization. The 
black curves show for comparison the general case stud- 
ied in Sees. IVTAl and lyiBl The dashed fines are for the 
same models but without tensors, i.e., ro = or e = 0. 

Let us first comment on Fig. 16 (amplitude 
parametrization). Since it was the correlation compo- 
nent that modified the acoustic peak structure in the 
generally correlated case, and this component is missing 
in the uncorrelated case, we do not have a mechanism 
that would favour a small ujc and large 6, Hq and Q\. 
Therefore these background parameters are now closer 
to their "adiabatic" values. Also smaller Uiso is now 
preferred. This is compensated by a slightly larger . 
Much larger ao is allowed but, as seen from the a-p plot, 
this leads to a smaller overall nonadiabatic contribution 
to Ci, since the correlation component is missing. The 
determination of ro is unaffected for the same reason as 
in the general case: isocurvature with a large niso modi- 
fies the small-scale end and the tensor contribution with 
nearly scale-invariant spectrum the large-scale end of the 



Ce spectrum. Taking the numbers from Table IV we find 
that, compared to the mixed model, the adiabatic model 
has odds 4.9 : 1 {Bf = -1-1.59) with tensors, and 3.9 : 1 
{Bf = +1.35) without tensors. Since the CMB data are 
compatible with the pure adiabatic model and the un- 
correlated mixed model has two parameters less than the 
generally correlated model, it is natural that the uncorre- 
lated model is much less disfavoured than the generally 
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Figure 16: Special cases, amplitude parametrization, 
CMB data, comparison of uncorrelated or fully cor- 
related models to the generally correlated model. The 

solid lines are for models with tensors. The black color is for 
the generally correlated model, the magenta dashed lines for 
fully anticorrelated (7 = —1), magenta dot-dashed for fully 
correlated (7 = +1), and blue for uncorrelated (7 — 0) model. 



correlated model that gave Bf +5 in support of the 
adiabatic model. 

Our result ai < 0.112 with tensors or 0.127 without 
tensors in Table |V] agrees well with the WMAP7 result 
[Tirro2] ai < 0.13 (their "Axion" model). For the nona- 
diabaticity we find ar < 3.7% with tensors and 4.3% 
without tensors. 



In the slow-roll parametrization (Figs. 17 and 18 1 e re- 



ceives a tighter constraint than i n th e general case, since 
the sin^ A factor is now 1 in Eq. (39 1. The constraint on 
To loosens. 

In the general case we had a tight connection between 
riiso and rias [see Eq. (39)] and the data (loosely) con- 
strained rtas- This led to a similar constraint on n-iso- 
However, now n^s is absent, and hence n-iso does not re- 

As the data seem to 
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ceive any constraint, see Fig. 
prefer the isocurvature contribution not showing up at 
any scale (which is achieved the easier the closer to 3 niso 
is; recall the left panel of Fig. 10 where riiso = 3.3), the 
posterior of n-isa bas a peak almost at the upper boundary 
of its prior range, niso = 1.15. 

The odds in favor of the corresponding adiabatic mod- 



Figure 17: Special cases, slow-roll parametrization, 
CMB data, comparison of uncorrelated or fully cor- 
related models to the generally correlated model. The 

line styles are the same as in Fig. |16[ 
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Figure 18: Special cases, slow-roll parametrization, 
CMB data, comparison of uncorrelated or fully cor- 
related models to the generally correlated model. The 

line styles are the same as in Fig. |16| 



els compared to the uncorrelated mixed model in the 
slow-roll parametrization are 17.6 : 1 {Bf — -f2.87) with 
tensors, and 80 : 1 {Bf = -f4.38) without tensors. 

Interestingly, the slow-roll parametrization gives very 
similar constraints on a — that are almost independent 
of the scale, since njgo ~ 1 — (ao < 0.110 with and 0.125 
without tensors) as the amplitude parametrization gives 
for ai. The 95% C.L. upper limits for ax are also very 
similar, although the shape of the posterior is different, 
the slow-roll parametrization preferring more clearly a 
zero nonadiabatic contribution. 
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B. Fully (anti)correlated perturbations 

In the fully correlated cases there is no Vnrik) spec- 
trum (7 = — 1 or 7 = +1). In the standard two-field 
inflation approach, we get this situation, when the curva- 
ture perturbations are completely dominated by the part 
that evolved from the entropy perturbations after horizon 
exit, and the tensor perturbations that were generated at 
horizon exit together with the uncorrelated part of the 
curvature perturbations are negligible. 

Therefore, in this case we assume no tensor perturba- 
tions, and we have just 8 independent parameters: back- 
ground parameters 



and 4 perturbation parameters 

In Al , rias , n-iso , "o ■ 



(67) 



(68) 



According to Eq. (39), the two-field slow-roll inflation 
gives 



= 1 + 277, 



(69) 



which reduces the number of independent perturbation 
parameters to 3 in the slow-roll parametrization. 

In the amplitude parametrization, the primary pertur- 
bation parameters are 



\nAl , lnyl2 I Oil : Cl2 : 

and in the slow-roll parametrization 



(70) 



(71) 



The marginalized 1-d posteriors for this model with the 
CMB data are indicated by the dashed magenta (7 = —1, 
100% anticorrelation) and dot-dashed magenta (7 = -1-1, 
100% correlation) curves in Fig. 16 for the amplidude 



parametrization, and in Figs. 17 and 18 for the slow-roll 
parametrization. 

As we would expect, in the amplitude parametriza- 
tion (see Fig. 16), the 7 = -f 1 case leads to very similar 



results as we got with the general correlation without 
tensors (dashed black curves). Some of the features of 
the well-fitting phenomenological mixed model are am- 
plified, such as the preference for very large niso. Due 
to the 100% correlation, for a given a the nonadiabatic 
contribution modifies the Ci spectrum in the maximal 
way. Therefore is formally constrained very tightly, 
but the peak of acoro moves to slightly larger value than 
in the general case. The posterior of ax is practically 
unaffected. On the other hand, the 7 = — 1 case differs 
significantly from the general case. Indeed, the ampli- 
tude parametrization results are now very close to the 
slow-roll parametrization results. This is natural, since 
the slow-roll parametrization favoured the negative cor- 
relation. The only clearly noticeable difference lies in 
the lower limit for ax that is now —0.08. Our result 



ai < 0.011 co-incides with the WMAP7 result [TIITB^. 
see their "Curvaton" model. 

The odds for the adiabatic model compared to the fully 
correlated mixed model are 206 : 1 {Bf — -1-5.33) and the 
fully anticorrelated model 11048 : 1 {Bj = +9.31). Here 
the overwhelmingly small odds of the mixed models are 
partially due to our choice to keep ai and 012 indepen- 
dent, i.e., allowing for riiso ^ fias- 

In the slow-roll parametrization this extra freedom is 
missing and we find that the adiabatic model has odds 
240 : 1 {Bf = +5.48) when compared to the fully corre- 
lated model, and 257 : 1 {Bf ~ +5.55) when compared 
to the fully anticorrelated mixed model. 



IX. COMPARISON TO OLDER RESULTS 



Some of the earliest constraints |194H197) on isocur- 
vature were obtained soon after COBE |198j in the mid 
1990s. Setting tight observational constraints on isocur- 
vature [199] became possible when Boomerang [200] and 
Maxima 2011 firmly observed the first acoustic peak in 
the CMB in the year 2000. Pre-WMAP constraints were 
derived, e.g., in [331 ISSl 12021 I203j . some papers focus- 
ing on the determination of the cosmological constant in 
the presence of primordial isocurvature modes jl87l I204| , 
while [205' finally ruled out a pure CDM isocurvature 
perturbation even in a spatially curved Universe. Con- 
straints from the Boomerang 2003 flight were derived in 
[206] and from Cosmic Background Imager (CBI) obser- 
vations in [207]. The WMAP first-year data [208J were 
used in [3 e [m m [Ml IMl [Wl mmTE\ to constrain 
various mixtures of adiabatic and isocurvature perturba- 



tions in a spatially flat Universe. The focus in |219j was 
in testing how much the (possible) presence of isocur- 
vature modes affects the determination of the geometry 
(spatial curvature) of the Universe. 

The WMAP 3-year data j220] were employed in ^ 
[T^ [Ti] and [T3] (note the post-publication corrections 
[httpTz/arxiv . org/pdf /astro-ph/ 0606685v3| 



to 



m 



for the case of a varying isocurvature spectral index). 
The most recent observational constraints and model se- 
lection with WMAP 5-year data combined with SN 
and SDSS matter power spectrum data were presented in 
[l6] for a spatially flat Universe. In [17] a comprehensive 
Bayesian comparison of flat and curved Universes both in 
pure adiabatic and in mixed cases was performed, finding 
that the spatial curvature is disfavoured roughly by the 
same amount as the isocurvature compared to the flat 
adiabatic model. 

In [222] the possibility of extra radiation (a radiation 
component other than photons or 3 standard species of 
neutrinos) carrying isocurvature is studied concluding 
that the current data allow the existence of an extra radi- 
ation component but do not favor its isocurvature mode. 
Ref. |223j shows that constraints on dark energy isocur- 
vature are very weak. 

Finally, forecasts for future experiments have been 
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worked out, e.g., in pTlE^lH^gUj . 

For a closer comparison to our results we pick six pub- 
lications where the pivot scale(s) and parametrizations 
are the most similar to our work. In our earlier work 
pS] . [5] and [T7], where we used 1-year, 3-year, and 
5-year WMAP data, respectively, and older small-scale 
CMB data, we did not include tensor perturbations. In 
the slow-roll parametrization this would correspond to 
assuming e = 0, but as all those work were in a phe- 
nomenological parametrization, we pick from Table [V] 
the "Ampl. par." results. Our new result (with CMB 
data) is ao < 0.158, while in [25] we obtained ao < 0.18; 
in [5 ao < 0.169; and in |T7j ao < 0.22. The actual 
limits on a are very sensitive to the precise amount of 
correlation preferred by the data, since the same a leads 
to a larger modification of the Ci spectrum if I7I dif- 
fers significantly from zero. The nonadiabatic contribu- 
tion to the CMB temperature variance is less sensitive 
to this "arbitrariness", and has evolved consistently to- 
ward tighter limits: from -0.075 < a^ < 0.075 pS] . 
0.017 < ar < 0.073 [9 , and -0.031 < ar < 0.067 [T?! 
to our new constraint —0.021 < ax < 0.053 without and 
-0.030 < ar < 0.049 with tensors. The WMAP 1-year 
data had relatively large error bars and did not extend 
to the third acoustic peak. As seen above, the range for 
ax was symmetric about 0. The later data map more 
precisely the acoustic peak structure and, as explained 
earlier, there is a persistent formal preference for a pos- 
itively correlated isocurvature component with a large 
riiso in the phenomenological approach. 

Of other publications, the most similar to our work are 
[T^ [T^ , although they do not perform Bayesian model 
comparison, but just parameter estimation. Ref. |14j 
uses WMAP 3-year data together with or without the 
SDSS DR4 LRG matter power spectrum, and includes 
tensor perturbations fixing the tensor spectral index nx 
by the first consistency relation. In addition, although 
their parametrization is phenomenological, they impose 
an "inflationary constraint" that riiso = ^i-as = 1 + nx- 
Therefore their work is closer to our slow-roll case than 
to our phenomenological approach: although we don't 



have riif 



in the slow-roll, our slow-roll approach 



guarantees that riiso is nearly scale invariant like theirs. 
Since they use spectral index parametrization with pivot 
scale k — 0.05 Mpc~^, we need to map our slow-roll re- 
sults to spectral indices and then from our pivot scale 
ko = 0.01 Mpc~^ to their pivot scale. Furthermore, we 
need to map to their isocurvature and correlation pa- 
rameters: Ba = \/a/(l — a) and cos 9a = sign(7)-y/|7[, 
respectively. After these mappings, our general slow- 
roll case with (without) tensors, using CMB data, gives 
Ba < 0.25 (0.26) at 95% C.L. and -0.79 < cosOa < 0.73 
(-0.67 < cosOa < 0.72). Adding the SDSS DR7 LRG 
data into the analysis, these change to Ba < 0.23 and 
-0.89 < cosOa < 0.21. Their resuhs with old WMAP 
and SDSS DR4 LRG are Ba < 0.28 (0.33) and -0.77 < 
cos < 0.27 (-0.40 < cos < 0.28). For tensor-to- 
scalar ratio at fc = 0.05 Mpc~^ we find r < 0.24 with 



CMB and 0.17 with CMB and MPK data, while their 
upper bound is 0.32. The general agreement is good 
when keeping in mind the differences in the data sets 
and parametrizations. For the special cases (uncorrelated 
model with tensors, uncorrelated model without tensors, 
fully anticorrelated, fully correlated) we find with CMB 
data Ba < 0.32, 0.35, 0.083, 0.077, while they find with 
the old WMAP data — in remarkable agreement — 0.31, 
0.33, 0.080, 0.087, respectively. For the tensor-to-scalar 
ratio at fc = 0.05 Mpc~^ in the uncorrelated case we find 
r < 0.25 while their constraint is 0.26. 

Ref. [E] uses the same CMB, SN and MPK data sets as 
we (except QUaD replaced by CBI and Boomerang) . Un- 
fortunately, they don't specify how the BAO part of the 
MPK theory and likelihood code was modified in order to 
take into account the isocurvature BAO when using the 
LRG sample, or whether the MPK code was not modified 
relying on the fact that BAO from subdominant isocurva- 
ture component would not change the results significantly 
with the current accuracy of data, as we have explicitly 
shown here. (BAO in the presence of isocurvature modes 
is further discussed in [24 ^ il69nl72j .) Ref. [19^ employs 
the phenomenological amplitude parametrization with- 
out tensor contribution and therefore the corresponding 
case in our study is "Ampl. par." (no tensors) with ei- 
ther CMB and SN or CMB and MPK data. As they have 
included both SN and MPK simultaneously, we would ex- 
pect their results to lie somewhere between our CMB&SN 
and CMB&MPK results, and be shghtly tighter. Indeed, 
taking an "average" of the dashed and dot-dashed curves 
in our Figs. [TT] and [121 noticing that athcir — 1 — ctour, and 
comparing to their figures 2 and 3 shows a good qualita- 
tive agreement, except for the exact shape of a2. Note 
that their cos A equals sign(7) They do not specify 

their sign convention, but we deduce from their figure 1 
that it is the same as ours. This is supported by their 
figure 4, (cos A, flm) plot, that indicates the same degen- 
eracy as we have seen, i.e., positive correlation prefers a 
large f^A- Finally, their numerical 95% C.L. constraint 
on ai (ao) is 6.1% (14.6%) while we find 6.7% (13.7%) 
with CMB and SN and 9.8% (16.7%) with CMB and 
MPK. We cannot do exact one-to-one comparison as they 
have one parameter less in their model, since they assume 

Finally, uncorrelated models are constrained in [18] 
which assumes — f^iso• Closest of our cases is the 
uncorrelated slow-roll model without tensors. However, 
even if njso is forced close to one in this model, it is 



completely free from riar, see Eq. (64 1 remembering that 
e = 0. As the CMB data are very "adiabatic" and hence 
require that the isocurvature contribution is subdomi- 
nant at all multipoles, they prefer spectral index riiso ^ 3. 
However, our slow-roll prior requires riiso ^ 1.15. There- 
fore the blue dashed curve in our Fig. [17] peaks at 
'liso = 1.15. This means that our model is quite dif- 
ferent from the model of [18] where the posterior peaks 
at riiso = n^T ~ 0.987. We find ao < 0.125. Map- 
ping into the spectral index parametrization and to pivot 
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scale k = 0.002 Mpc"! used in [TH] we obtain a < 0.117, 
whereas their constraint 0.19 with QUaD and WMAP 
5-years data is much looser, showing — once again — 
that the formal constraints on undetected isocurvature 
parameters can depend strongly on the theoretical as- 
sumptions, here whether one imposes riar — n\so or uses 
slow-roll parametrization that makes them independent. 



X. CONCLUSIONS 

Apart from using more recent CMB data, the main 
new results of this paper compared to |17j . which also 
utilizes MultiNest and Bayesian model comparison, are: 
1) Showing that allowing for a primordial tensor contri- 
bution has a negligible effect on the determination of the 
nonadiabatic contribution and that the possible presence 
of a CDM or baryon isocurvature mode has a minor ef- 
fect on the determination of the tensor-to-scalar ratio, 
as long as the tensor spectral index obeys the first infla- 
tionary consistency relation. 2) Showing that the results 
change considerably compared to jTTj if, instead of the 
phenomenological parametrization, one uses inflationary 
slow-roll parametrization that forces all the spectra near 
to scale invariance. 3) Including one of the most impor- 
tant complementary (to CMB) data sets into the anal- 
ysis, namely the matter power spectrum, and showing 
that its current accuracy is not enough to give direct ex- 
tra constraints on isocurvature. However, even the cur- 
rent MPK data change the results considerably due to 
the tension between the low matter density preferred by 
the CMB data and the high matter density preferred by 
the MPK data. 4) In [H] (and [H]) all the runs were 
done without CMB lensing, because already computa- 
tionally intensive runs turn by an order of magnitude 
heavier if lensing is applied. Ignoring the lensing was 
justified by the fact that the Id posteriors stayed rather 
unaffected. However, in this paper, in Sec. |VD] we have 
noticed that comparing the unlensed mixed model to the 
unlensed adiabatic model may make the mixed model 
seem slightly more favourable than it actually is, in par- 
ticular what comes to the Bayesian evidences. Therefore, 
the CMB lensing is now included in the analysis. In the 
end the effect on model comparison is mild: for exam- 
ple, Bf = +5.48 with and Bf — -1-4.18 without lensing 
when comparing the adiabatic model to the mixed one 
(phenomenological model with tensors). 

We conclude that current cosmological data give no 
evidence of an isocurvature perturbation contribution. 
Thus we can only set upper limits to it. Since isocur- 
vature perturbations contribute to the large-scale CMB 
more strongly than to small-scale CMB, in comparison to 
the adiabatic perturbations, the tightest limits to isocur- 
vature perturbations are at the large scales. In the ab- 
sence of a detection of an isocurvature contribution, the 
phenomenological approach to constraining the parame- 
ters then leads to a preference of a steeply blue primordial 
isocurvature perturbation spectrum (large Uiso) unless a 



prior constraint on spectral indices is imposed. The limits 
to an isocurvature contribution then come mainly from 
the relative locations of the acoustic peaks in the CMB 
angular power spectrum. 

In addition to a phenomenological approach, we con- 
sidered an approach motivated by two-field inflation, 
where the four first-order slow-roll parameters were fit- 
ted to the data, with the prior assumption that the slow- 
roll parameters are small. This means that all primor- 
dial spectra are constrained to be close to scale invari- 
ant. Since in this case isocurvature perturbations con- 
tribute much more to large CMB scales than to small 
scales, when compared to the adiabatic perturbations, 
the isocurvature contribution is mainly constrained by 
the relative levels of the low-£ and high-£ parts of the 
observed CMB angular power spectrum. This leads to 
a tighter limit to isocurvature perturbations than when 
their spectral index is not constrained a priori. When 
spectra are close to scale invariant, the contribution of 
isocurvature perturbations to Ci are somewhat similar 
to that of tensor perturbations — both contribute mainly 
to low £. If isocurvature perturbations are anticorrelated 
with adiabatic perturbations, they have a negative contri- 
bution at low £, and thus their contribution can be partly 
canceled by tensor perturbations. Thus the combination 
of negatively correlated isocurvature perturbations and 
tensor perturbations fit the data better than positively 
correlated isocurvature perturbations and tensor pertur- 
bations. However, since for a given isocurvature fraction 
a, correlated perturbations give a much stronger signal 
in the CMB than uncorrelated perturbations, allowing 
for tensor perturbations does not relax the upper limit 
to a. 

In the two-field slow-roll infiation approach it is not 
feasible to go beyond the level of first order in slow-roll 
parameters, since the number of parameters would then 
become too large. In any case, the slow-roll approxima- 
tion is only valid, when the parameters are small. There- 
fore we assumed a prior constraint that their magnitudes 
are smaller than 0.075. With this prior assumption, we 
obtained a somewhat tighter posterior limit to the slow- 
roll parameter related to the first derivatives of the infia- 
tion potential V{(j)i, (j)2), namely e < 0.048, but no useful 
constraint to the second-derivative slow-roll parameters 
rjij. Thus to make progress on constraining two-field in- 
flation, we need more accurate data, e.g., from the Planck 
satellite. 
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